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ABSTRACT 


Fluid  shear  stress  is  one  of  the  most  important  quantities  charac¬ 
terizing  the  interaction  between  a  flowing  fluid  and  a  solid  body.  A 
knowledge  of  surface  shear  stress  is  essential  in  determining  the  drag 
on  submerged  bodies,  the  energy  losses  in  internal  flows,  local  surface 
heat  transfer  rates  in  both  internal  and  external  flows,  flow  separation, 
cavitation,  and  unsteady  flow  behavior  in  rotating  machinery. 

In  spite  of  many  previous  Investigations  of  both  laser  Doppler 
anemometry  (LDA)  and  surface  hot-film  gauges,  the  use  of  these  methods 
for  determining  wall  shear  stress  is  poorly  developed.  Conventional 
LDA  measurements  near  surfaces  are  plagued  by  spurious  background  scattering 
and  marginal  resolution.  Use  of  surface  hot-film  gauges  is  severely  in¬ 
hibited  by  the  necessity  of  empirical  calibration  and  an  incomplete  under¬ 
standing  of  their  response  in  turbulent  and  unsteady  flows. 

This  research  program  is  developing  new  methodology  for  the  use  of 
thes^  -tyvo  complementary  techniques  as  quantitative  measurement  tools  for 
determ+hing  local  surface  shear  stress.  The  LDA  system  is  based  on 
detecting  laser-induced  fluorescence  from  solid  particles  in  a  so-called 
"backscatter"  mode.  Two  "'/idual  velocity  components  can  be  measured 
by  using  multiplexing  in  ’equency  domain,  the  hot-film  probes  employ 

multiple  eteni&nts  such  that  ventional  calibration  procedures  can  be 

eliminated.  Both  theory  and  experiment  are  used  to  determine  the  utility 
of  LDA  and  hot-film  instrumentation  in  unsteady  flows. 
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I.  DESCRIPTION  OF  PROGRAM 

A.  Rationale 

The  vast  majority  of  all  energy  conversion  processes  involve  fluid 
flow  and  heat  transfer.  The  performance  of  gas  turbines,  compressors, 
heat  exchangers,  and  other  equipment  used  to  generate  power  cannot  be 
optimized  without  extensive  theoretical  and  experimental  research.  And 
a  key  element  in  this  research  is  the  ability  to  make  accurate  measurements 
of  detailed  flow  characteristics  within  the  devices  being  tested. 

The  importance  to  the  ONR  Power  Program  of  fluid  flow  and  hear  trans¬ 
fer  research,  including  research  on  measurement  techniques,  may  be  judged 
by  the  numerous  research  projects  being  supported  in  these  areas.  A  series 
of  workshops  have  been  sponsored  by  Project  SOUID  to  assess  progress  in  par¬ 
ticular  aspects  of  fluid  flow  research;  since  1969,  three  conferences  have 
explored  gas  dynamics  and  unsteady  flow  phenomena,  two  have  been  devoted  to 
turbulence,  and  three  have  had  the  theme  of  laser  instrumentation  for  fluid 
flow  measurement.  As  one  reads  the  Proceedings  of  these  workshops  [e.g. 

Murthy  (1976)],  it  becomes  clear  that  complicated  three-dimensional  unsteady 
flows  are  of  utmost  importance;  and  our  ability  to  make  accurate  measurements 
in  such  environments  is  severely  limited. 

Fluid  shear  stress  is  one  of  the  most  important  experimental  quantitites 
characterizing  the  interaction  between  a  flowing  fluid  and  a  solid  body.  A 
knowledge  of  surface  shear  stress  is  essential  in  determining  the  drag  on  sub¬ 
merged  bodies,  the  energy  losses  in  internal  flows,  local  surface  heat  trans¬ 
fer  rates  in  both  external  and  internal  flows,  flow  separation,  cavitation, 
and  unsteady  flow  behavior  in  rotating  devices.  High-frequency  fluctuations 
of  shear  stress  can  pinpoint  acoustic  modes  in  inlets  and  ducts,  or  detect 
unsteady  shock-induced  boundary  layer  separation.  Local  shear  stress  (both 
magnitude  and  direction)  is  a  quantity  that,  in  some  cases,  may  be  predicted 
by  theory;  and  comparison  between  theory  and  experiment  for  this  quantity  is 
a  sensitive  test  of  the  accuracy  of  the  calculation  procedure. 
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There  are  only  two  practical  methods  of  measuring  the  fluid  shear 
stress  acting  on  a  variety  of  solid  surfaces:  laser  Doppler  anemometry 
(LDA)  to  determine  the  fluid  velocity  profile  near  the  wall,  thereby 
obtaining  the  shear  rate  and,  with  knowledge  of  the  fluid  viscosity, 
the  shear  stress;  and  the  surface  hot-film  gauge  whose  rate  of  heat 
loss  can  be  related  to  local  shear  stress  by  appropriate  calibration. 

Other  methods  (e.g.  Preston  Tube,  hot-wire  anemometers,  electrochemical 
probes,  and  floating-element  balances)  have  inherent  limitations  that 
preempt  their  use  in  most  practical  circumstances. 

B.  Objectives 

The  objectives  of  our  research  program  are  as  follows: 

(1)  Develop  a  comprehensive  first-principles  description  of  the 
thermal  response  of  surface  hot-film  gauges  to  imposed  fluid  flow.  The 
theory  will  include  complete  3-dimensional  transient  conduction  to  the 
substrate  and  a  time-dependent  quadratic  fluid  velocity  profile.  Experi¬ 
mental  verification  will  be  obtained  for  quartz-backed  and  poly  (vinyli- 
dine  fluoride)  lithographed  gauges  in  our  pulsatile  flow  facility  and  cone- 
plate  calibration  apparatus. 

(2)  Develop  a  two-component  backscatter  LDA  system  with  fluorescent 
seed  particles  for  optimum  S/N  ratio  and  improved  spatial  resolution  near 
solid  surfaces.  Bragg-shift  multiplexing  of  a  single  excitation  wavelength 
will  be  used  to  provide  two  orthogonal  velocity  components.  Make  a  direct 
comparison  between  2-component  LDA  using  fluorescence  and  conventional  2- 
component  system  using  two  different  wavelengths. 

(3)  Produce  a  detailed  quantitative  comparison  between  LDA  and  hot- 
film  instruments  in  standard  laminar  and  turbulent  flow  regimes  (pulsatile 
flow,  pipe  flow,  cone-plate  laminar  and  turbulent  flow).  Refine  interpre¬ 
tation  of  the  hot-film  response  on  the  basis  of  measured  turbulent  sub-layer 
velocity  components  including  intermittent  sweeps. 

The  shear  stress  measurement  techniques  that  we  are  developing  will 
have  substantial  potential  value  to  a  wide  variety  of  problems  of  interest 
to  the  Navy.  Measurements  within  complex  heat  exchanger  geometries  will 


provide  valuable  data  for  understanding  local  heat  transfer  coefficients. 

Measurements  can  be  made  on  turbine  and  compressor  models  to  detect 
flow  separation  and  to  unravel  complex  flow  velocity  patterns  that  defy 
theoretical  or  conventional  experimental  methods.  Measurements  may  also 
be  easily  made  on  large  ship  models  and  full-scale  vessels  of  local  fric¬ 
tion  drag. 

(4)  Fabricate  multiple-film  gauges  using  microlithography  and  thick- 
film  technology  to  achieve  a  self-calibrated  sensor  with  sufficient  infor¬ 
mation  to  give  shear  direction  as  well  as  magnitude. 

(5)  Apply  these  composite  hot-film  gauges  to  high-speed  unsteady 

flow. 

(6)  Develop  quantitative  assessments  and  recommendations  regarding 
relative  ease  of  use,  accuracy,  and  cost  of  LDA  vs.  hot-film  technology 
for  measuring  wall  shear  stress  in  a  variety  of  flows  of  interest. 

C.  Personnel 

The  following  personnel  participated  in  the  program  during  the  current 
contract  period: 

C.  Forbes  Dewey,  Jr.  -  Co-Principal  Investigator 

Peter  W.  Huber  -  Co-Principal  Investigator 

Kenneth  M.  Kalumuck  -  Research  Assistant 

Don  F.  Pardo  -  Research  Assistant 

Leanna  Tyau  -  Research  Assistant 

Judi  Steciak  -  Research  Assistant 

Richard  Struve  -  Research  Assistant 

Marj  Joss  -  Secretary 

R.  Fenner  -  Research  Staff  Member 


D.  Major  Accomplishments  for  Current  Reporting  Period 

Fig.  1  is  a  table  summarizing  important  previous  theoretical  calculations 
of  hot-film  probe  heat  transfer  to  a  flowing  fluid.  The  current  theory  being 
developed  by  Ken  Kalumuck  is  substantially  more  complete  in  representing  the 
conditions  existing  in  practice  than  any  previous  efforts  in  this  area.  The 
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work  of  Pardo,  completed  with  partial  support  from  the  current  program, 
in  our  laboratory,  was  Important  In  establishing  the  quantitative  impor¬ 
tance  of  substrate  conduction. 

Kalumuck  considers  a  thin  film  on  a  wall  surface  (the  x-z  plane) 
exposed  to  a  flowing  fluid.  The  film  is  of  arbitrary  shape  and  an  amount 
of  heat  q(x,z)  per  unit  area  is  generated  by  the  film.  Some  heat  is  con¬ 
ducted  into  the  wall,  and  additional  heat  is  given  up  directly  to  the 
fluid. 

A  solution  to  the  coupled  heat  conduction  and  fluid  convection 
equations  is  obtained.  First,  the  equations  are  subjected  to  Fourier 
transformation  in  two  of  the  three  spatial  directions.  This  leads  to  two 
ordinary  differential  equations  in  y  (the  coordinate  normal  to  the  wall). 

The  solution  to  these  equations  is  then  obtained  in  terms  of  the  transformed 
variables.  As  a  final  step,  a  two-dimensional  Fast  Fourier  Transform  (FFT) 
is  performed  to  obtain  the  detailed  temperature  field  in  any  desired  x-z 
plane.  The  resulting  solution  includes  all  of  the  detailed  information  on 
heat  flux  both  to  the  solid  and  to  the  fluid.  Specific  results  have  been 
obtained  and  are  reported  in  Section  II. 

Preliminary  results  have  been  obtained  by  R.  Struve  using  an  ion 
laser  to  excite  fluorescent  particles  and  using  the  fluorescence  to 
track  the  particles.  From  these  experimental  results  and  the  calculations 
of  L.  Tyau,  we  expect  that  our  technique  of  using  fluorescent  particles  in 
the  LDA  experiments  will  be  successful.  The  details  of  these  results  will 
be  found  in  Section  III. 


T, 
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II.  HOT-FILM  GAUGES  FOR  MEASURING  FLUID  SHEAR  STRESS 


A.  Review  of  Literature 

Hot-film  elements  [Ludweig  (1950)]  have  been  used  for  many  years  to 
measure  flow  velocities  and  wall  shear  in  a  variety  of  fluid  flow  appli¬ 
cations.  Methodology  important  to  their  use  was  first  developed  by 
Lieomann  and  Skinner  (1955)  and  Ling  (1955).  Important  contributions  to 
the  theoretical  interpretation  of  the  relation  between  heat  and  momentum 
transfer  for  such  probes  was  made  by  Lighthill  (1954),  and  Ling  (1963), 
and  subsequently  by  Brown  (1967),  Spence  and  Brown  (1968),  Ackerberg  et  al . 
(1978),  Springer  and  Pedley  (1973),  and  (using  a  mass  transfer  analogy) 

Hanratty  and  colleagues  [Fortuna  and  Hanratty  (1962,  1963),  and  Sirkar  and 
Hanratty  (1970a,  1970b)].  A  current  review  by  Hanratty  and  Campbell  is 
also  in  press  [Hanratty  (1982)]. 

Among  the  more  important  examples  of  the  use  of  hot-film  techniques 
are  those  of  Ling  et  al_.  (1968,  1969),  Bellhouse  and  Schultz  (1967,  1968), 

Brech  and  Bellhouse  (1973),  Clark  (1974),  Armistead  and  Keyes  (1968), 

Eckelmann  and  Reichardt  (1971),  Pope  (1971),  Geremia  (1970),  McCroskey 
and  Durbin  (1972),  and  Murthy  and  Rose  (1978). 

The  principle  of  operation  of  a  hot-film  gauge  is  simple:  the  rate 
of  removal  of  heat  from  a  small  heated  film  is  related  to  the  instantaneous 
shear  stress  (or  more  precisely,  the  velocity  gradient)  in  the  fluid  immediately 
adjacent  to  the  film.  Important  experimental  difficulties  and  theoretical 
limitations  to  the  accurate  quantitative  use  of  such  gauges  have  been  iden¬ 
tified  previously  and  are  discussed  in  the  following  sections.  But  the 
potential  advantages  of  hot-film  gauges  over  existing  alternatives  are  clear. 
Hot-film  gauges  are  physically  small  and  require  no  more  wiring  or  Deripheral 
instrumentation  close  to  the  measuring  point  than  does  a  pressure  transducer 
or  conventional  strain  gauge.  They  can  operate  in  opaque  liquids  or  at 
optically  inaccessible  points  within  rotating  machinery,  piping  systems  and 
valves.  They  can  be  designed  to  withstand  corrosive  or  otherwise  hostile 
environments.  Flush  mounting  on  a  wall  can  assure  that  the  measurement 
process  will  result  in  minimal  disturbance  to  the  fluid  phenomena  under 


investigation.  The  gauges  are  potentially  inexpensive  and  simple  to  use; 
that  they  cannot  he  described  in  this  way  at  present  can  be  attributed  pri¬ 
marily  to:  (1)  difficulties  encountered  in  accurate  calibrations;  and  (2) 
substantial  uncertainty  with  regard  to  the  interpretation  of  the  hot-film 
signal  in  unsteady  and  turbulent  flows.  The  promise  of  this  instrument  is 
nonetheless  clear;  the  hot-film  gauge  could  become  as  versatile  and  easy  to 
use  in  measuring  rates  of  strain  (or  shear  stresses)  in  fluids  as  are  con¬ 
ventional  strain  gauges  for  measuring  ordinary  strains  in  solids. 

The  underlying  Drinciple  of  operation  of  thin  film  shear  stress  gauges 
stems  from  the  Reynolds  analogy  between  heat  transfer  and  local  skin  fric¬ 
tion  coefficient  in  laminar  two-dimensional  boundary  layers  in  which  heating 
effects  due  to  friction  and  fluid  compression  are  negligible.  In  its  gen¬ 
eral  form  the  Reynolds  analogy  emerges  directly  from  the  boundary  layer 
momentum  and  energy  equations  [Schlichting  (1968)]  and  states  that  the  heat 
transfer  to  or  from  a  wall  is  related  to  the  shear  stress  at  the  wall  by  a 
function  of  the  form: 

Nu  =  0.5  cf  •  Re  •  *  ,  Pr)  (1) 


where  $  represents  a  functional  dependence  that  includes  the  upstream 
history  of  the  boundary  layer.  Here  the  Nusselt,  Reynolds  and  Prandtl 
numbers  and  the  skin  friction  coefficient  are  defined  in  the  usual  way: 


Nu  = 


q* 


«VT-) 
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Re  = 
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In  Eqs.  (1)  -(5),  x  and  y  are  coordinates  of  position  parallel  to  and 
perpendicular  to  the  wall,  q  is  the  local  heat  transfer  rate  per  unit 
area,  i  is  a  characteristic  length  scale  in  the  flow  direction,  k  is 
the  fluid  thermal  conductivity,  Tw  is  the  local  wall  temperature,  Uw,  p^, 
and  Tm  are  the  free  stream  fluid  velocity,  density  and  temperature,  and 
v  and  a  the  fluid  kinematic  viscosity  and  thermal  diffusivity  (k/pCp) 
respectively;  (3u/3y)w  denotes  the  local  velocity  gradient  at  the  wall. 

The  theory  for  the  operation  of  thin  film  gauges  follows  directly 
from  Eq.  (1).  If  we  consider  the  case  of  a  thin  heated  film  of  length  L 
and  width  W  mounted  flush  with  an  adiabatic  wall,  with  the  film  tempera^ 
ture  Tw  (Figure  1),  then  for  any  given  laminar  flow  velocity  profile 


Figure  1.  Laminar  flow  over  a  heated  strip. 


above  the  film  it  is  possible  to  correlate  the  rate  of  heat  transfer 
from  the  film  to  the  fluid  with  the  shear  stress  at  the  wall.  For  the 
case  of  a  velocity  profile  that  can  be  considered  linear  over  a  distance 


from  the  wall  larger  than  the  thermal  boundary  layer  thickness  above 
the  probe,  the  functional  relation  in  Eq.  (1)  can  be  determined  exactly 
[Leveque  (1928)].  Averaging  the  heat  transfer  over  the  film  surface  for 
this  case  one  obtains 


The  average  heat  transfer  rate  per  unit  area,  q,  may  be  equated  to  the 
electrical  dissipation  in  the  probe  divided  by  the  size  of  the  probe, 
i  .e. 


(9) 


where  V  is  the  voltage  across  the  thin  film,  R  is  the  film  resistance,  and 

L  and  M  are,  respectively,  the  film  dimensions  in  the  direction  of  flow  and 

normal  to  it.  It  is  frequently  assumed  that  the  two-dimensional  analysis 

leading  to  Eq.  (5)  is  valid  if  L  «  W;  this  is  not,  in  fact,  the  case. 

In  practice,  Eq.  (6)  does  not  hold  because  the  wall  is  not  adiabatic 

★ 

and  heat  conduction  occurs  into  the  solid.  This  leads  to  heating  of  the 
fluid  in  the  boundary  layer  upstream  of  the  probe  and,  in  general,  a  trans¬ 
fer  of  heat  from  the  heated  fluid  into  the  wall  downstream  of  the  probe. 


For  moderately-conducting  solids,  the  thermal  penetration  depth  into  the 
solid  is  comparable  to  the  probe  length  W,  and  the  heat  conduction  Drob- 
lem  is  3-dimensional . 
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Eq.  6  suggests  that  the  relation  between  the  measured  electrical  heat 
dissipation  and  shear  stress  is  of  the  form 

%  ■ A  * B (w)'  m 

where  A  and  B  are  constants  to  be  found  by  calibration,  NuM  is  the  measured 
Nusselt  number  obtained  by  substituting  Eq.  (9)  into  Eq.  (2),  and  AT  is 
assumed  to  be  fixed  for  both  the  calibration  and  the  conditions  of  use. 

Before  proceeding,  it  should  be  recognized  that  the  flow-independent 
term,  "A”,  in  Eq.  (10)  is  often  2  to  4  times  larger  than  the  term  con¬ 
taining  information  on  the  local  shear  stress.  This  term  arises  primarily 
from  substrate  conduction  [Tanner  (1967)]  and  is  relatively  insensitive 
to  natural  convection  [Ackerberg  et al_.  (1978)].  Measurement  sensitivity 
is  increased  by  decreasing  substrate  conduction  (reducing  "A"  relative  to 
"B  (3u/9y)w"h  but  this  imposes  certain  difficulties  with  regard  to  the 
interpretation  of  measurements  in  unsteady  flows. 

Our  experience  [c.f.  Sdougos  (1977)]  as  well  as  the  experience  of 
most  other  investigators  [c.f.  Pope  (1971),  Bellhouse  and  Schultz  (1966), 
Geremia  (1970),  Brown  and  Davey  (1971),  and  Murthy  and  Rose  (1978)]  is 
that  Eq.  (10)  represents  an  excellent  fit  to  laminar  steady-flow  cali¬ 
bration  experiments  under  conditions  where  the  parameter  NL  is  moderately 
large,  say  >  20.  For  smaller  values  of  N^,  the  theory  of  Ackerberg 
et  at.  (1978)  may  be  consulted.  In  flows  with  pressure  gradients,  for 
which  the  region  of  linear  velocity  variation  is  reduced  (Fig.  2),  addi¬ 
tional  terms  in  Eq.  (10)  are  necessary  [see  Spence  and  Brown  (1968)]. 

The  fact  remains,  however,  that  conventional  use  of  hot-film  shear 
stress  gauges,  even  in  steady  laminar  flow,  requires  calibration  under 
conditions  closely  approximating  those  of  the  proposed  measurements. 

This  is  an  extremely  tedious  and  cumbersome  procedure.  One  principal 
objective  of  our  research  is  to  substantially  reduce  (effectively  elimin¬ 
ate)  the  tedious  calibration  procedures  that  are  presently  required. 


There  is  very  little  quantitative  information  available  in  the 
literature  on  the  response  of  hot-film  gauges  to  unsteady  shear. 
Although  the  electrical  response  time  of  the  film  and  feedback  circuit 
is  many  KHz,  the  gauge  cannot  respond  in  a  quasi -steady  manner  to  the 


Figure  2.  The  thermal  boundary  layer  growing 
into  a  linear  shear  flow. 


imposed  shear  stress  fluctuations  unless  the  convective  heat  transfer  and 
the  substrate  conduction  heat  transfer  approach  a  steady  state  within  a 
time  comparable  to  that  between  fluctuations. 

To  understand  the  nature  of  this  problem,  we  consider  the  response  of 
a  heated  film  to  flow  fluctuations  imposed  upon  it  by  the  external  flow. 

The  easiest  case  to  consider  is  a  flat  plate,  upon  which  a  thin  film  gauge 
is  mounted,  that  oscillates  parallel  to  its  surface.  In  a  coordinate  sys¬ 
tem  moving  with  the  surface,  this  is  equivalent  to  an  oscillating  boundary 
layer  immediately  above  the  plate  —  a  problem  treated  by  Lighthill  (1954), 
Bellhouse  and  Schultz  (1968),  and  Fortuna  and  Hanratty  (1971).  Experimental 
evidence  is  available  from  Fortuna  (1970),  Ling  et_  al_.  (1968)  and,  for  free 
probes,  Clark  (1974).  We  wish  to  present  a  simple  physical  argument  for 


estimating  the  upper  frequency  limit  for  which  the  steady-state  relation 
of  Eq.  (1)  between  wall  shear  stress  and  measured  heat  transfer  is  valid. 

It  will  then  be  shown  that  this  simple  model  reproduces  the  essential  fea¬ 
tures  of  the  more  rigorous  solution  of  Fortuna  and  Hanratty  (1971). 

Referring  to  Fig.  2,  we  define  u^  to  be  the  fluid  velocity  at  a  height 
6y  above  the  surface.  It  is  clear  that  the  thermal  boundary  layer  developing 
along  the  strip  will  not  be  quasi-steady  unless  the  fluid  particles  heated 
by  the  strip  have  nearly  constant  velocity  during  their  passage  past  the 
strip.  This  suggests  that  a  particle  moving  along  a  streamline  6y  above 
the  plate  should  travel  a  distance  of  at  least  10L  during  one  cycle  of 
the  imposed  velocity  fluctuation  of  frequency  in  order  to  insure  quasi¬ 
steady  behavior;  i.e. 

Wp  =  2frf  <  2ir  ^  j  .  (11) 


The  magnitude  of  u^  for  a  probe  operating  in  the  linear  velocity  regime 
(Fig.  2)  is  equal  to  the  velocity  gradient,  (3u/3y)w  ,  times  the  height, 
6y  ,  of  the  thermal  layer  at  the  rear  of  the  strip.  From  Ling  (1963), 


~  2.9 

L  "  Nl!/3 


(12) 


where 
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L20u/3y)w 
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Pr  L*2 


(13) 


L*  =  Lu*/v 

[v(3u/3y)w  ]l/2 


(14) 


u*  = 


Inserting  Eqs.  (12)  and  (13)  into  (11),  we  have  the  restriction  on  fre¬ 
quency  for  quasi -steady  probe  response: 


This  result  may  be  cast  into  a  form  suitable  for  comparison  with 
the  results  of  Fortuna  and  Hanratty  (1971).  We  define  a  nondimensional 
frequency 


By  using  Eqs.  (13)  -(14),  we  have  the  simple  form  for  the  present  model: 


1.8 


(17) 


Fortuna  and  hanratty  have  calculated  the  ratio.  A2,  of  the  mean  square 
fluctuation  amplitude  of  the  heat  transfer  to  a  flow  disturbance  of  fre¬ 
quency  u)  to  the  amplitude  of  the  induced  heat  transfer  which  would  be 
measured  if  the  flow  disturbance  had  the  same  amplitude  but  the  oscilla¬ 
tion  was  very  slow.  Their  result  is: 
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For  the  criteria  we  presented  In  Eq.  (17),  <*>£  L*2^3  <  1.8,  their  calcu¬ 
lations  predict  about  a  20%  decrease  In  the  mean  square  fluctuation,  or 
a  10%  decrease  In  the  (linear)  amplitude  of  the  heat  transfer  fluctuation. 
Referring  to  Eq.  (6),  the  implied  wall  shear  stress  fluctuation  would  then 
be  too  low  by  30%. 

Some  typical  numbers  are  informative  here.  For  a  pipe  flow  at  a 
characteristic  velocity  V  =  0.1  ms-1  in  a  pipe  of  0.01  m  diameter  with 
water  as  the  liquid  we  have: 
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For  a  probe  of  length  L  =  1.5  x  10" 4  m,  L*  -  1.35  and  from  Eq.  (14)  the 
upper  limit  for  the  frequency  response  is  found  to  be:  f  2  (<d/2tt)  <  10  Hz. 
This  typical  calculation  points  out  the  substantial  limitation  to  quasi¬ 
steady  frequency  response  that  hot  film  gauges  possess. 

The  quasi-steady  response  criteria  presented  in  Eq.  (15)  was  obtained 
by  considering  the  heat  transfer  from  the  film  to  the  flow.  Inasmuch  as 
a  substantial  portion  of  the  heat  generated  in  the  film  is  conducted  into 
the  substrate,  quasi-steady  behavior  also  requires  that  a  time-independent 
substrate  temperature  distribution  has  been  achieved.  The  time  scale 
for  this  equilibration  can  be  estimated  from  well-known  solutions  to 
unsteady  heat  conduction  problems  in  solids  [c.f.  Carslaw  and  Jaeger 
(1947)].  The  characteristic  penetration  depth,  xchar,  of  an  oscillating 
thermal  wave  into  a  solid  is  given  by 


X  /  2 

xchar  * 

where  tc$  =  (k/pc)  is  the  thermal  diffusivity  of  the  solid.  If  we  asso¬ 
ciate  x(.|iar  with  the  probe  length  L  and  the  time  t  with  a  half-period 
of  the  thermal  oscillation  [u>s  *  (ir/t)]  in  the  substrate,  then  the  sub¬ 
strate  temperature  distribution  will  be  in  equilibrium  provided 

“s  <  ^  lVL^  O9) 

The  rati*  siting  frequencies  calculated  from  fluid  response  and  sub 
strata  is  then 


(20) 


One  of  the  objectives  of  our  program  is  to  obtain  accurate  relations 
for  the  frequency  response  of  hot-film  probes.  Although  the  arguments 
presented  in  the  preceding  paragraphs  establish  limits  for  the  quasi-steady 
frequency  response  of  the  gauges  to  imposed  shear  stress  fluctuations, 
it  is  not  necessary  to  limit  their  use  to  such  low  frequencies  provided 
that  the  gauge  transfer  function  is  known  as  a  function  of  frequency. 

As  Eq.  (19)  demonstrates,  this  transfer  function  must  include  the  effects 
of  substrate  conduction. 

Previous  attempts  to  establish  such  transfer  functions  by  empirical 
measurements  [Armlstead  and  Keyes  (1968a,  1968b),  Bellhouse  and  Schultz 
(1966),  Seed  and  Wood  (1970),  Woodruffe  (1970),  Clark  (1974)]  have  been 
inconclusive  at  best.  In  particular,  Clark  (1974),  Smith  et  at.  (1974) 
and  Bellhouse  and  Schult2  (1966,  1967)  observed  a  transfer  function  that 
exceeds  unity  at  intermediate  frequencies,  and  rolls  off  approximately  as 
uf1  at  high  frequencies.  This  result  is  in  contradiction  to  the  monotoni- 
cally  decreasing  transfer  function  predicted  by  Fortuna  and  Hanratty  (1971). 
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Clark  (1974)  has  remarked  that  the  theory  of  Pedley  (1972)  [see  also 
Pedley  (1976)]  might  be  applicable,  but  this  theory  does  not  predict  the 
correct  phase  relationship  between  heat  transfer  and  shear  stress  and  is 
therefore  suspect. 

We  anticipate  that  our  research  will  encourage  use  of  hot-film 
gauges  in  unsteady  flows  where  a  mean  bulk  or  core  flow  has  a  periodic  or 
fluctuating  component  characterized  by  a  single  frequency.  If  a  boundary 
layer  of  thickness  6m  due  to  the  mean  flow  separates  the  wall  from  the 
core  flow  region,  periodic  fluctuations  in  the  core  flow  will  be  felt  as 
shear  stress  fluctuations  at  the  wall  only  if  they  are  below  some  maximum 
frequency.  This  result  is  the  classical  problem  of  the  "depth  of  penetra¬ 
tion"  of  a  viscous  wave  [Stuart  (1955)]:  the  viscous  boundary  layer  is 
impervious  to  shear  stress  fluctuations  at  frequencies  higher  than  or  of 
order: 


(21) 

m 

For  a  boundary  layer  thickness  10”*  m  and  a  kinematic  viscosity  of  10~6 
m2  s"1  the  frequency  cut-off  is  of  order  100  Hz.  More  will  be  said  about 

the  implications  of  ■''his  result  in  the  paragraphs  that  follow. 

The  simplest,  and  at  this  stage  probably  the  only,  possible  argument 
for  assessing  the  performance  of  hot-film  gauges  in  turbulent  flows  is  that 
when  the  thermal  boundary  layer  thickness  is  less  than  that  of  the  laminar 
sublayer  the  principles  and  equations  outlined  and  developed  in  the  previous 
section  should  be  applicable.  The  well-known  linear  velocity  predictions 
for  a  sublayer  extending  a  distance  of  order  5  v/u*  above  the  wall  are 
usually  cited,  and  the  requirement  that  the  thermal  boundary  layer  be  thinner 
than  this  is  thus  of  the  form  [Ling  (1963);  Liepmann  and  Skinner  (1954)]: 


(22) 
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In  terms  of  the  non-dimensional  variables  usually  used  to  describe  the 
turbulent  boundary  layer, 

U+  =  (u/uj  (23) 

y+  s  (yu*/v)  (24) 

\  =  [v(3u/8y)Jl/2  (8) 

this  criteria  becomes 

S j  <  5  (25) 

Estimating  6^  using  Eq.  (12),  we  have  a  criterion  for  the  maximum  length, 

L,  of  a  hot-film  gauge  that  would  be  suitable  for  measurements  of  shear 
stress  in  turbulent  flow: 

L*  =  (Lu*/v)  <5.1  Pr  (26) 

This  criterion  is  normally  satisfied  with  the  small  values  of  L  (10“**  m  or 
less)  normally  employed  in  hot-film  gauges.  Such  arguments,  as  far  as 
they  go,  cannot  be  faulted;  but  there  is,  as  yet,  no  proof  that  the 
Reynolds  analogy  holds  in  the  laminar  sublayers  of  turbulent  flow  when 
Eq.  (26)  Is  satisfied.  It  seems  safe  to  expect  that  under  some  condi¬ 
tions  turbulent  flow  wall  shear  stresses  inferred  from  hot  film  data 
based  on  the  conventional  laminar  boundary  Tayer  equations  will  be  rea¬ 
sonably  accurate.  On  the  other  hand,  enough  is  already  known  about  the 
three  dimensional  structures  and  periodic  "sweeping"  of  laminar  sublayers 
to  suggest  that  parallels  between  ordinary  laminar  boundary  layers  and 
laminar  sublayers  are  not  perfect.  An  objective  of  our  research  program 
is  to  define  the  relationship  between  mean  and  fluctuating  shear  stresses, 
as  measured  with  a  hot-film  gauge,  and  the  true  dynamic  processes  occurring 
in  the  sublayer  of  a  turbulent  boundary. 
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The  important  measurements  that  have  been  made  with  hot-film  gauges 
in  turbulent  flow  will  now  be  reviewed. 

Bellhouse  and  Schultz  (1968)  describe  the  measurement  of  laminar 
oscillatory  shear  flows  superposed  on  large  mean  flows  in  air  with  a  flush- 
mounted  hot-film  gauge  (constant  temperature  mode).  The  authors  suggested 
that  the  probe,  when  operating  in  the  laminar  sublayer  of  a  turbulent  flow, 
would  respond  as  it  would  in  a  truly  laminar  rectilinearly  oscillatory 
flow,  but  they  did  not  offer  any  evidence  to  support  this  claim. 

Armistead  and  Keyes  (1968a,  1968b)  compare  the  frequency  spectra  of 
heat  transfer  from  a  flush-mounted  hot-film  in  a  turbulent  flow  to  the 
spectra  of  wall  mass  transfer  and  wall  pressure  fluctuations.  The  authors 
expect  that  fluctuations  in  the  rates  of  heat  and  mass  transfer  and  wall 
pressure  fluctuations  "should  be  closely  related"  since  all  are  induced 
by  "flow  perturbations  in  the  wall  layer."  This  is  not  correct;  Townsend 
(1976)  remarks  that  wall  pressure  fluctuations  are  influenced  by  velocity 
fluctuations  within  a  distance  (larger  than  the  sublayer  thickness)  com¬ 
parable  to  the  wavelength  of  the  velocity  fluctuations.  In  fact,  Armistead 
and  Keyes  do  note  that  the  intensity  of  pressure  fluctuations  is  roughly 
two  orders  of  magnitude  larger  than  heat  and  mass  transfer  intensities. 

The  overall  comparison  between  the  three  measures  of  wall  fluctuation 
(heat  and  mass  transfer  and  wall  pressure)  given  by  Armistead  and  Keyes 
(1968a,  Fig.  8)  does  not  support  their  claim  of  closely  related  fluctua¬ 
tions  of  these  three  measures.  Also,  the  authors'  claim  for  the  applica¬ 
bility  of  wall  hot-film  probes  to  detect  dynamic  wall  shear  stress  using 
the  "demonstrated  relationship  between  mean  rate  of  heat  transfer  and  the 
wall  shear  stress"  is  not  verified  by  their  data. 

Experiments  by  Bakewell  and  Lumley  (1967)  were  conducted  to  produce 
large  sublayer  thicknesses  within  which  data  on  streamwise  velocity 
fluctuations  could  be  obtained.  From  the  correlations  obtained,  a  flow 
pattern  was  proposed  which  consisted  of  randomly  distributed  counter¬ 
rotating  eddy  pairs  of  elongated  streamwise  extent  (Fig.  3).  The  authors 
conclude  that  this  nonlinear  vortex  stretching  mechanism  makes  invalid 
earlier  suggested  theories  of  linearized  sublayer  behavior.  The  implica- 
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tions  of  this  theory  with  regard  to  hot-film  response  has,  to  date,  not 
been  addressed. 

The  earlier  work  of  Mitchell  and  Hanratty  (1966)  and  the  theoretical 
study  by  Fortuna  and  Hanratty  (1971)  overlook  the  sublayer  picture  sug¬ 
gested  by  the  model  of  Bakewell  and  Lumley.  The  Hanratty  group  uses  mass 
transfer  techniques  operating  with  very  large  Schmidt  numbers  (~5000, 


Figure  3.  Sublayer  model  of  Bakewell  and  Lumley  (1967). 

typically)  and  it  was  supposed  that  the  only  significant  sublayer  motions 
affecting  mass  transfer  In  the  sublayer  were  the  longitudinal,  recti¬ 
linear  oscillations:  they  explicitly  neglected  velocities  normal  to  the 
flow  direction  and  they  also  assumed  instantaneous  longitudinal  velocities 
were  linear  in  distance  from  the  wall. 

Later  work  by  Sirkar  and  Hanratty  (1970a,b)  acknowledged  the  impor¬ 
tance  of  the  sublayer  picture  emerging  from  Bakewell  and  Lumley' s  studies 
since  the  contribution  to  time-dependent  mass  transfer  due  to  longitudinal 
flow  fluctuations  is  found  to  be  negligible  in  comparison  to  that  due  to 
velocity  gradients  in  the  cross-stream  direction  and  the  direction  normal 
to  the  wall.*  In  this  event,  the  mass  transfer  probe  cannot  be  used  to 
measure  wall  shear  stress  fluctuations.  Kline  et  ai.  (1967)  provide  fur¬ 
ther  evidence  that  the  flow  in  the  sublayer  is  spatially  organized,  time- 
dependent  large  sublayer  eddies  providing  definite  regions  of  updraft  and 
sidedraft,  as  in  Lumley's  model,  with  rather  small  flow  fluctuations  in 
the  longitudinal  direction. 

~  --r  -  >  —  ' 

This  conclusion  has  yet  to  be  verified! 
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One  question  that  has  been  addressed  by  previous  investigators  is 
the  following:  Is  the  calibrated  relation  between  hot-film  heat  transfer 
and  wall  shear  stress  obtained  in  laminar  flow  applicable  to  a  turbulent 
boundary  layer?  As  Pope  (1971)  points  out,  the  original  data  of  Llepmann 
and  Skinner  (1954)  were  not  accurate  enough  to  resolve  this  question. 

Pope  (1971)  derives  a  correction  based  on  local  turbulence  intensity  that 
appears  to  make  the  laminar  and  turbulent  calibrations  coincide;  but  very 
limited  data  are  available  to  test  the  universality  of  this  correction, 
and  the  theoretical  basis  of  the  correction  appears  to  be  highly  question¬ 
able.  Brown  (1967)  and  Bellhouse  and  Schultz  (1966)  find  good  agreement 
between  laminar  and  turbulent  flow,  but  the  data  are  very  limited.  In 
supersonic  external  airflow,  where  laminar  sublayers  are  very  thick  in 
comparison  with  internal  flows,  Owen  and  Bellhouse  (1970)  found  good 
agreement  between  laminar  and  turbulent  calibrations  for  the  mean  shear 
stress  as  a  function  of  gauge  heat  loss.  No  evidence  exists  regarding 
the  effects  of  temperature  fluctuations  on  either  gauge  calibration  or 
the  measurement  of  fluctuating  quantities. 

Only  two  examples  are  known  to  us  of  the  use  of  wall  gauges  to  mea¬ 
sure  the  fluctuating  component  of  shear  stress.  Fortuna  and  Hanratty 
(1971),  using  mass-transfer  gauges  analogous  to  heat-transfer  gauges, 
measured  the  r.m.s.  fluctuation  in  wall  shear,  to  be  0.35  t^.  Their 
measurements  were  corrected  for  frequency  response  as  discussed  in  a 
previous  paragraph.  Eckelmann  (1974)  used  a  wall-mounted  hot-film  gauge 
and  measured  r'  *  0.24  t  .  This  agreed  very  closely  with  the  result  he 

W  W 

obtained  by  appropriate  interpretation  of  hot-wire  velocity  measurements. 
Assuming  that  the  instantaneous  velocity  profile  near  the  wall  is  linear, 
then  the  instantaneous  velocity  u  and  shear  stress  t  are  related  to  the 
instantaneous  fluctuations  if  and  ?' 


u  =  u  +  if 

(27) 

and 

r  -  7  +  r 

(28) 
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He  measured  the  rms  velocity  fluctuation  u‘  to  be  0.24  u  for  values  of 
y  £0.1  (i.e.  well  within  the  conventional  laminar  sublayer  boundary  of 
y  =5).  Therefore,  his  measurements  of  u'  yield  r'  =  0.24  7,  in  agree¬ 
ment  with  the  hot-film  result.  By  an  independent  method,  Popovich  and 
Hummel  [(1967a,  1967b);  see  also  Popovich  (1969)]  obtained  x'  =  0.30  7W< 

A  further  independent  estimate*  may  be  made  by  combining  wall  pressure 
fluctuation  data  summarized  by  Leehey  (1966)  and  reported  by  Blake  (1970) 
with  skin  friction  correlations  measured  by  Burton  (1974).  The  result 

of  this  calculation  is  t*  =  0.22  7  . 

w  w 

The  measurements  of  Eckelmann  were  performed  under  most  favorable 
conditions;  the  laminar  sublayer  of  his  oil  channel  was  1  cm  thick,  10-100 
times  thicker  than  is  available  with  most  facilities  and  sufficiently  large 
that  his  hot-film  wall  gauge  had  a  length  L*  [see  Eq.  (13)]  of  only  about 
0.03.  We  expect  that  the  ~  50%  difference  between  the  results  of  Fortuna 
and  Han  ratty  (1971)  and  Eckelmann  (1974)  arise  in  part  from  the  supra - 
1  inear  increase  of  u‘  with  distance  from  the  wall. 

Our  research  in  this  area  should  establish  a  much  firmer  foundation 
for  the  use  of  hot-film  gauges  to  measure  fluctuations  in  wall  shear  stress 
under  a  variety  of  experimental  conditions.  The  fact  remains  +hj.t  under 
normal  circumstances  the  hot-film  gauge  is  subject  to  cross-stream  shear 
fluctuations  as  well  as  streamwise  shear  fluctuations,  a  the  nature  of 
the  gauge  response  may  depend  upon  the  corns?  ition  between  these  two  quan¬ 
tities,  as  suggested  by  tne  Bakewell  and  Lumley  (1967)  model.  The  two 
central  questions  are:  1)  Do  three-dimensional  effects  vitiate  the  two- 
dimensional  boundary  layer  assumption  that  underlies  the  Reynolds  analogy?; 
2)  Do  transient  disruptions  or  distortions  of  the  laminar  sublayer  so  com¬ 
plicate  the  interpretation  of  the  hot-film  gauge  performance  that  its 
effective  use  in  turbulent  applications  is  precluded?  Relevant  to  both  of 
these  questions  is  the  feasibility  of  developing  very  simple  calibration 
schemes  that  properly  account  for  such  effects. 

In  the  past  decade,  a  significant  fraction  of  the  fundamental 
research  on  turbulent  boundary  layers  has  been  devoted  to  understanding  this 


We  are  Indebted  to  Prof.  P.  Leehey  of  M.I.T.  for  suggesting  this  calculation. 
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bursting  phenomenon  and  integrating  these  conceptual  models  with  existing 
data  on  average  and  statistical  properties  of  turbulent  flows.  Well  over 
100  papers  exist  on  this  subject,  and  it  is  beyond  the  scope  of  the  present 
discussion  to  treat  them  in  detail.  It  suffices  to  note  that  this  is  a 
difficult  problem  and  much  research  will  continue  in  this  area  for  many 
years  to  come.  It  is  clear,  however,  that  both  the  hot-film  gauge  and 
the  laser  Doppler  anemometer  are  potentially  very  powerful  tools  in  the 
continuing  investigation  of  this  problem. 

B.  Theoretical  Formulation 

The  model  we  are  developing  is  designed  to  deal  with  probes  of 
arbitrary  geometry,  and  to  account  for  convection  in  the  fluid  and  diffu¬ 
sion  in  all  directions  in  both  the  fluid  and  the  substrate.  The  probe  is 
modeled  as  a  heat  flux  source  of  arbitrary  (but  specified)  heat  flux 
distribution  and  geometry  located  on  a  plane  conducting  wall  (Fig.  4). 

The  fluid  above  the  probe  flows  with  a  uniform  shear  (i.e.,  a  linear  velo- 
* 

city  profile  ).  The  governing  equations  are  the  three-dimensional  fluid  and 
solid  energy  equations: 

v  =  l29) 

V*2T$*  =  0,  y*  <  0  (30) 


It  will  be  possible  to  incorporate  various  other  velocity  profiles  into 
the  model.  The  linear  profile  is,  however,  an  excellent  assumption  for 
flows  in  which  the  thermal  boundary  layer  is  much  thinner  than  the  momen¬ 
tum  boundary  layer.  Such  will  be  the  case  for  most  liquids  and  for  flows 
in  which  the  thermal  boundary  layer  begins  its  growth  downstream  of  the 
beginning  of  the  momentum  boundary  layer. 
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where 

v*2  =  _ii  +  _ii_  +  _iL  , 

3x*2  3y*2  3z*2 

s  =  fluid  velocity  gradient  (u*  =  sy*), 
u  =  fluid  velocity  in  the  x  direction, 

T  *  temperature  measured  above  ambient, 
a  =  fluid  thermal  diffusivity, 

*:  denotes  dimensional  quantities, 

f:  as  a  subscript  denotes  fluid  quantities,  and 

s:  as  a  subscript  denotes  solid  (substrate)  quantities. 

The  temperature  field  is  continuous  at  the  fluid-solid  interface.  The 
heat  flux  is  also  continuous  at  this  interface,  except  at  the  probe 
which  acts  as  a  source  of  heat.  Thus: 

T  *  =  T  *  at  y*  =  0  (31) 

f  s  J 

3T*  lr>\ 

k  ~-L-  -kf  — ~  =  Q*  (x*,  z*)  at  y*  =  0  w 

3y*  3y* 

where  k  is  the  thermal  conductivity  and  Q*  is  the  arbitrary  (but  specified) 
probe  heat  flux  distribution.  At  sufficiently  large  distances  from  the 
probe,  the  temperature  must  attain  the  ambient  value,  or 

Tf*,  T$*  -*■  0  as  x*.  y*  or  z*  ±  »  (33) 
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These  equations  can  be  cast  in  a  form  more  amenable  to  analysis 
by  appropriate  nondimensionalization.  We  define 

L  =  /a/s  a  characteristic  length  for  thermal  diffusion 

=  the  total  power  output  of  the  probe  (=  J  J  Q*dx*dz*) 

and  introduce  dimensionless  quantities: 

x  =  x*/L,  y  =  y*/L,  z  =  z*/L 

T  -  T*  kfL/Ps 

K=ks/kf 

Q ( x »  z)  =  Q*(x*,  z*)  •  L2/Ps. 

Equations  (29) -(33)  can  now  be  expressed  as; 


v2Tf  -  y 

d'f. 

3x 

=  0,  y  >  0 

(34) 

V2Ts  =  0 

»  y 

<  0 

(35) 

Tf  =  V 

y 

=•  0 

(36) 

3Ts 

K  ay 

3Tf 

3y 

=  Q(x,  z),  y  =  0 

(37) 

Tf,  is  0 

as  x,  y,  z  ±  » 

(38) 

A. 
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The  solution  of  Eqs.  (34)  -  (38)  depends  on  only  the  conductivity 

ratio  K  and  the  source  distribution  Q(x,  z).  The  shear  s  is  part  of  the 

definition  of  a  characteristic  diffusion  length  (L  =  /a/s ) .  A  flow  Peclet 

number  can  be  obtained  from  the  definition  of  Q:  for  a  characteristic 

z 

streamwise  probe  dimension  of  length  L*  ,  the  Peclet  number  is  SL*  /a  = 

,  2  P  P 

U*/L)  . 

Equations  (34)  -  (38)  represent  a  general  model  for  analyzing  hot- 
film  probes  in  steady  uniform  shear  flows.  We  are  developing  an  efficient 
and  flexible  solution  procedure  for  this  set  of  equations  which  retains  all 
terms  of  the  equations,  allows  for  all  ranges  of  the  values  of  the  governing 
parameters,  and  permits  the  specification  of  arbitrary  probe  geometries  and 
heat  fluxes  (i.e.,  arbitrary  Q(x,  z),  including  the  analysis  of  multiple- 
film  probes.  Solution  of  this  system  of  equations  yields  a  relation  between 
probe  heat  flux  and  the  resulting  temperature  field.  A  numerical  "calibration" 
of  a  probe  design  is  then  made  possible  by  knowledge  of  the  relations  between 
fluid  shear,  probe  temperature,  and  probe  heat  flux. 


Method  of  Solution 

The  temperature  field  at  y  =  0,  Tq  (x,  z),  is  obtained  from  our  model 
[Eqs.  (34)  -  (38)]  by  a  combination  of  analytical  and  numerical  techniques. 
Solving  for  the  temperature  field  at  only  one  value  of  y  (y  =  0,  here) 
greatly  increases  computational  efficiency.  Solutions  to  the  energy  equations 
in  both  the  fluid  [Eq.  (34)]  and  solid  [Eq.  (35)]  domains  are  obtained  by 
first  analytically  performing  Fourier  transforms  in  both  the  streamwise  and 
spanwise  directions.  These  two-dimensional  transformations  reduce  the  par¬ 
tial  differential  equations  (34^and  (35)  to  ordinary  differential  equations 
whose  solutions  are  then  found.  The  arbitrary  constants  of  these  solutions 
are  evaluated  from  the  Fourier  transforms  of  the  boundary  conditions  [Eqs. 

(36)  -  (38)].  This  method  of  solution  is  similar  to  that  used  by  Tanner  [1967] 
and  by  Davies  and  Kimber  [1972]  in  solving  simplified  forms  of  the  governing 
equations. 

The  two-dimensional  Fourier  transform  of  the  wall  temperature  dis¬ 
tribution  is  defined  as 


T0  (p,  q) 


00 


e2ir1(px  *  T0  (x,  2)  dxdz. 


•00 


00 


(39) 
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We  have  found  the  solution  of  Eqs.  (34)  -  (38)  in  two-dimensional  transform 
space  to  be: 


ri/6 


Ai(fl)/Ai '(n) 


T.  ( p *  q)  ~  0  tP»  q)  ■  ]/o  i/p 

0  (2*p)  '  Ai(Q)/Ai’(^)  -  1 

n  =  e*i/3  (2tt)4/3  fiLi 

P‘ 7 


(40) 


where 


A.(fi)  =  the  Airy  function  of  complex  argument  ft, 
A1'(a)  =  (Ai(fl)), 

p  =  transform  variable  in  the  x  direction, 
q  =  transform  variable  in  the  z  direction,  and 

A 

Q  =  the  two-dimensional  transform  of  Q. 


The  desired  temperature  field,  Tq  (x,  z),  is  then  found  by  application 
of  the  two-dimensional  inverse  Fourier  transform: 


T0  (x.  2)  = 


(px  *  ,2)  T  (p>  q)  dpdq  . 


(41) 


All  of  the  information  specifying  the  details  of  the  probe  geometry 


and  its  heat  flux  distribution  is  contained  in  Q  (p,  q).  The  analyses 
of  most  other  researchers  can  be  viewed  as  special  cases  of  this 
solution.  For  example,  an  adiabatic  wall  solution  is  obtained  by  setting 
the  conductivity  ratio  K  to  zero.  A  two-dimensional  geometry  can  be 
analyzed  by  making  the  heat  flux  distribution  Q  (x,  z)  very  much  longer 
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in  the  spanwise  direction  than  in  the  streamwise  direction. 

The  complexity  of  the  integrand  (40)  of  inversion  integral  (41) 
necessitates  a  numerical  evaluation.  The  inversion  is  achieved  by 
applying  a  series  of  numerical  one-dimensional  Fast  Fourier  Transforms 
(FFT's)  first  in  one  direction  and  then  in  the  other.  The  integrand  (40) 
to  which  the  FFT's  are  applied  must  first  be  evaluated  numerically  at  a 
finite  set  of  discrete  values  of  (p,  q).  The  Airy  function  and  its  deri¬ 
vative  are  evaluated  from  series  expansions.  Taylor  series  are  used  for 
small  values  of  |n|  ,  and  asymptotic  series  are  used  for  larger  values  of 

A 

the  argument.  The  transform  of  the  heat  flux  source  distribution,  Q  (p,  q), 
can  be  found  analytically  for  relatively  simple  source  distributions 
(e.g.,  uniform  flux  over  a  rectangular  probe),  and  can  be  obtained  numeri¬ 
cally  for  more  complicated  distributions.  All  ooerations  are  performed 
using  complex  arithmetic. 

The  numerical  evaluation  of  Eqs.  (40)  and  (41)  is  currently  under 
development.  The  initial  cases  being  investigated  include  a  point  source 
of  heat  mounted  on  an  adiabatic  wall  and  a  square  uniform  heat  flux  probe 
exposed  to  a  moderate  Peclet  number  flow  and  mounted  on  a  wall  of  variable 
conductivity.  Preliminary  calculations  of  the  wall  temoerature  distributions 
along  x  =  0  and  along  z  =  0  for  a  square  uniform  heat  flux  probe  mounted  on 
a  adiabatic  wall  at  a  Peclet  number  of  36  are  shown  in  Fig.  5.  These  cal¬ 
culated  distributions  are  qualitatively  reasonable:  along  z  =  0  the  peak 
temperature  occurs  above  the  probe  but  downstream  of  its  midpoint,  the 
temperature  distributions  decay  monotonically  in  all  directions  from  the 
maximum  value,  and  the  streamwise  temperature  profiles  are  skewed  such  that 
the  downstream  decay  is  much  slower  than  the  upstream  decay.  The  spanwise 
temperature  profiles,  as  expected,  are  symmetric  about  z  =  0. 

Streamwise  (z  =  0)  temperature  profiles  for  a  square  uniform  heat 
flux  probe  mounted  on  an  adiabatic  wall  at  Peclet  numbers  ranging  from 
4  to  64  are  presented  in  Fig.  6.  Here,  the  temperature  profiles  are  nor¬ 
malized  by  multiplying  by  the  probe  length,  L  .  The  point  of  maximum 
temperature  moves  further  upstream  with  increasing  Peclet  number.  As  the 
Peclet  number  increases,  the  decay  of  the  temperature  distributions 
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proceeds  less  rapidly  downstream  and  more  rapidly  upstream.  The  average 
probe  temperature  is  also  seen  to  decrease  with  increasing  Peclet  number 
(indicating  a  probe  heat  transfer  rate  per  unit  average  temperature  that 
increases  with  Peclet  number).  The  effect  of  varying  conductivity  ratio 
on  the  streamwise  (z  =  0)  temperature  profiles  due  to  a  square  uniform 
heat  flux  probe  at  constant  Peclet  number  (here  Peclet  number  equals  16) 
is  presented  in  Fig.  7.  Increasing  the  conductivity  ratio  results  in 
calculated  temperature  profiles  that  are  lower  in  magnitude,  whose  maxi¬ 
mum  move  closer  to  the  probe  center,  and  that  are  more  symmetric  about 
these  maxi  mums. 

That  the  calculated  temperature  distributions  have  not  decayed  to 
zero  results  from  the  use  of  a  discrete  rather  than  continuous  transform. 

This  results  in  a  truncation  of  the  physical  space  computational  domain. 
Another  possible  source  of  error  in  the  computed  results  is  the  use  of  a 
finite  rather  than  infinite  transform.  Tests  are  being  conducted  to  de¬ 
termine  the  discretization  and  truncation  which  will  introduce  negligible 
error  in  the  computed  temperature  distributions.  An  appropriate  average 
integrated  value  of  the  integrand  at  (p,  q)  =  (0,  0)  has  been  found  to 
deal  with  the  singularity  there. 

C.  Continuing  Research 

Testing  and  refinement  of  our  computational  model  is  currently  in 
progress.  We  are  generating  solutions  for  simple  probe  geometries  and 
flow  situations  to  be  checked  against  existing  analytical  and  numerical 
solutions  and  available  experimental  data.  After  this  initial  verifica¬ 
tion  and  refinement,  we  will  investigate  more  comolicated  probe  and  flow 
geometries  with  the  dual  goal  of  gaining  a  better  understanding  o^'  the 
complicated  heat  transfer  processes  involved  and  of  conducting  parametric 
P^obe  design  studies.  Such  analyses  will  enable  us  to  provide  a  "numerical 
calibration"  cf  a  probe  design.  Experimental  work  utilizing  hot-film  probes 
for  shear  stress  measurement  in  both  steady  and  unsteady  flows  has  been  in 
progress  in  our  laboratory  for  some  time.  A  future  stage  of  our  experimental 
effort  will  be  to  fabricate  probes  of  different  geometries  and  substrate 


Figure  7.  Effect  of  Varyinq  Conductivity  Rat 
Profiles  for  a  Square  Probe;  Pe=16 


conductivities  which  have  been  designed  by  our  computational  model. 

These  can  then  be  tested  under  a  variety  of  flow  conditions  enabling 
us  to  verify  cur  "numerical  calibrations"  and  to  gain  additional 
understanding  of  the  underlying  physics.  Our  steady  flow  model  and 
experiments  should  also  provide  much  needed  insight  into  the  performance 
of  probes  in  unsteady  flows. 
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III.  LASER  DOPPLER  ANEMOMETRY 


A.  Review  of  Literature 


The  laser  Doppler  anemometer  (LDA)  has  undergone  extensive  develop¬ 
ment  in  the  past  decade  and  is  now  a  versatile  and  accurate  instrument  for 
measuring  velocities  with  minimal  disturbance  to  the  fluid  flow  under  in¬ 
vestigation.  Several  excellent  monographs  on  the  theory  and  operation  of 
the  LDA  are  available  [Durst,  Mel  ling  and  Whitelaw  (1976),  Durrani  and 
Greated  (1977)];  and  a  recent  review  paper  [Buchave,  George  and  Lumley 
(1979)]  describes  the  state  of  the  LDA  art  in  the  measurement  of  tur¬ 
bulence  . 

The  LDA  optical  configuration  that  is  now  used  most  frequently  is 
termed  the  "differential  Doppler"  or  "fringe"  mode  of  operation;  it  is 
illustrated  in  Figure  8.  Two  intersecting  laser  beams  form  an  interference 


Figure  8. 
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fringe  pattern  in  the  fluid  that  is  orthogonal  to  the  component  of 
velocity  to  be  measured.  In  ths  absence  of  a  Bragg  cell  frequency 
shifter,  this  fringe  pattern  is  stationary.  Any  small  particle  that 
passes  through  the  fringe  pattern  will  scatter  light,  producing  a  sig¬ 
nal  at  the  detector  with  a  frequency  proportional  to  the  particle  velo¬ 
city  normal  to  the  fringes  and  inversely  proportional  to  the  fringe 
spacing. 

If  the  fringes  are  stationary,  the  detector  will  not  distinguish 


between  particles  moving  in  opposite  directions  with  the  same  velocity. 

By  introducing  a  Bragg  cell  in  one  of  the  optical  paths,  the  beams  can 
be  made  to  have  slightly  different  frequencies  and  the  resultant  fringe 
pattern  will  be  a  moving  one,,  the  fringes  continuously  sweeping  across 
the  probe  volume  in  one  direction.  The  scattered  light  signal  from  a 
particle  moving  through  the  probe  will  now  be  sensitive  to  the  direction 
of  the  particle  motion.  A  further  refinement  is  the  use  of  multi-colored 
or  polarized  beams  and  additional  detectors  to  measure,  simultaneously, 
two  or  three  components  of  the  particle  velocity. 

It  is  not  our  intention  to  present  here  a  comprehensive  review  of 
the  experimental  challenges  that  arise  in  the  many  and  varied  applications 
of  LDA  instrumentation.  We  propose  to  study,  compare,  and  calibrate  against 
each  other,  two  methods  of  measuring  wall  shear  stress  in  unsteady  and  tur¬ 
bulent  flows;  and  we  discuss  here  only  those  factors  of  LDA  operation  that 
are  most  important  in  this  context. 

The  way  in  which  the  LDA  can  be  used  to  determine  wall  shear  stress 
is  conceptually  straight-forward.  The  two  orthogonal  components  of  fluid 
velocity  parallel  to  the  surface  of  a  wall  are  measured  with  the  LDA  at  a 
series  of  locations  above  the  surface.  If  (x,  z)  are  the  coordinates  in 
the  plane  of  the  surface  and  (u,  w)  are  the  velocities  in  these  directions, 
then  the  magnitude  of  the  surface  vector  shear  stress  for  a  Newtonian  fluid 
is 

i1/2 

,  (42) 

where  y  is  measured  normal  to  the  surface  and  (3u/3y)w  and  (9w/3y)w  are 


/  9u  \  2  , 

/  aw  \  , 

(w  ay/wi 

r  > 


-40- 


the  measured  velocity  gradients  at  the  wall. 

The  accurate  determination  of  ?w  and  the  fluctuating  shear  stress 
t'  require  a  careful  assessment  of  potential  errors  and  the  use  of  appro- 

W 

priate  corrections  to  the  raw  measurements.  These  corrections,  summarized 
by  Buchhave  et  al.  (1979),  include  the  following: 

(1)  A  correction  for  so'-rious  frequency  components 
arising  from  the  entry  and  exit  of  individual  particles 
from  the  measuring  volume.  This  term  is  referred  to  as 
"transit  time  broadening";  it  is  not  present  when  the  so- 
called  "burst  mode"  or  single-particle  analysis  is  used. 

(2)  A  correction  for  bias  arising  from  the  presence 
of  velocity  gradients  across  the  measuring  volume  (gra¬ 
dient  broadening"). 

(3)  A  correction  to  account  for  the  phase  fluctua¬ 
tions  arising  from  correlation  of  the  light  scattered 
from  several  particles  within  the  measuring  volume 
["turbulence  broadening";  see  George  and  Lumley  (1973)]. 

In  addition  to  these  calculable  and  correctable  sources  of  broadening,  the 
experimenter  is  faced  with  the  necessity  of  achieving  good  spatial  resolu¬ 
tion,  designing  a  stable  and  repeatable  optical  system  that  can  traverse 
a  region  of  high  shear  with  an  accuracy  of  a  few  urn,  and  having  a  collec¬ 
tion  system  for  the  scattered  radiation  that  exhibits  a  very  high  signal- 
to-noise  ratio.  The  system  we  propose  to  use  has  these  attributes.  To 
understand  the  advantages  of  the  system  we  propose,  a  review  of  previous 
experiments  is  necessary. 

Table  1  presents  the  conditions  of  measurement  of  the  seven  previous 
published  experiments  we  consider  most  relevant  to  our  proposed  research. 

Of  these,  the  measurements  of  Karpuk  and  Tiederman  (1976)  and  Quigley  and 
Tiedeman  (1977)  are  the  most  successful  in  achieving  good  spatial  resolution 
near  the  wall  of  a  turbulent  flow.  These  excellent  experiments  were  able 
to  provide  confirmation  of  certain  of  the  correction  factors  predicted  by 
theory  [see  the  summary  of  Buchhave  et  al.  (1979)]  that  were  described  In 
the  previous  paragraph.  The  authors  demonstrate  that,  even  under  favorable 
circumstances,  uncorrected  LDA  signals  will  yield  apparent  mean  velocities 
10-30%  greater  than  the  true  mean;  and  measured  rms  fluctuations  may  be  in 
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Table  1.  LOA  Measurements  Near  Solid  Surfaces 
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* 

error  by  100%  or  more  .  After  appropriate  correction,  they  obtain  values 

of  7  that  are  in  agreement  with  pressure  drop  measurements  and  values  of 
"  + 

(x'/T  '  at  y  =  2  that  are  in  substantial  agreement  with  the  careful  mea- 
w  w 

surements  of  Eckelmann  (1974)  that  were  cited  in  a  previous  section. 

The  experiments  of  Reich  (1973)  and  Oldengarm  et  al.  (1975)  are 
somewhat  qualitative,  inasmuch  as  Reich  fails  to  account  for  significant 
transit  time  broadening  and  Oldengarm  et  al.  make  no  attempt  to  correct 
their  measurements  for  bias  arising  from  gradient  broadening.  Agrawal 
et  al.  (1978)  were  concerned  with  skewed  laminar  flow  in  -urves  tubes, 
and  no  attempt  was  made  to  make  detailed  measurements  near  the  tube  walls; 
furthermore,  their  laminar  flow  had  much  less  severe  velocity  gradients 
than  generally  occur  in  turbulent  or  highly  unsteady  flows.  As  far  as 
the  experiments  of  Yoganathan  et  al.  (1979)  are  concerned,  they  made  no 
effort  to  quantify  the  errors  of  their  measurements;  we  calculate  that 
their  reported  wall  shear  stresses  are  potentially  in  error  by  a  substantial 
amount. 

We  now  turn  to  the  question  of  signal -to-noise  ratio  of  the  scat¬ 
tered  radiation  received  by  the  detector.  A  careful  examination  of  various 
alternative  optical  arrangements  has  convinced  us  that  the  backscatter 
arrangement,  as  used  by  Born  et  al.  (1978)  and  others,  is  by  far  the  most 
useful  for  measurements  near  solid  surfaces.  The  angle  6  between  the  two 
intersecting  beams  must  be  large  to  achieve  small  probe  volume  dimensions; 
and  this  allows  relatively  large  collection  angles  for  the  backscattered 
radiation,  thereby  increasing  the  amount  of  received  signal. 

One  formidable  difficulty  in  all  IDA  measurements  near  surfaces  is 
that  imperfections  in  the  surface  itself  will  cause  unwanted  scattered 
radiation  to  be  collected  by  the  detector.  This  background  signal  is 
frequently  very  large  compared  to  the  amount  of  radiation  scattered  by  a 


This  •’S  consistent  with  the  large  errors  in  mean  velocity  observed  by 
Born  et  al.  (1978)  near  the  wall  of  a  small  tube,  although  the  correc¬ 
tions  Born  et  al.  compute  are  in  error.  Their  Eq.  (2)  contains  several 
misprints  [compare,  e.g.  Berman  and  Dunning  (1973)].  Further,  their 
computed  correction  does  not  agree  with  the  nroper  relation  of  Berman 
and  Dinning. 
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moving  particle  in  the  fluid.  An  example  of  how  this  background  scatter¬ 
ing  can  dominate  the  detector  signal  is  given  by  the  results  of  Oldengarm 
et  al.  (1975).  Normally  high-pass  electronic  filters  are  used  to  reduce 
the  magnitude  of  this  background,  but  such  electronic  discrimination  is 
marginal  at  best  and  cannot  remove  the  photocurrent  shot  noise  that  con¬ 
tributes  to  the  filtered  signal. 

An  elegant  method  of  removing  these  large  background  signals  was 
proposed  by  Stevenson  et  al.  (1975).  They  used  test  particles  containing 
fluorescent  dye.  As  the  test  particles  traverse  the  measuring  volume, 
the  incident  laser  light  is  absorbed  and  reradiated  as  fluorescent  emis¬ 
sion  at  longer  wavelengths.  Inasmuch  as  the  fluorescent  lifetime  of 
suitable  dyes  is  on  the  order  of  1-10  nsec,  the  fluorescent  intensity  is 
directly  proportional  to  the  instantaneous  fringe  field  intensity  exper¬ 
ienced  by  the  measuring  particle. 

Laser-induced  fluorescence  therefore  shifts  the  signal  radiation 
to  longer  wavelengths,  allowing  the  use  of  simple  optical  filters  to  remove 
all  unwanted  scattered  light  at  the  excitation  wavelength.  In  the  back- 
scatter  mode,  the  fluorescence  signal  is  comparable  to  that  which  would  be 
produced  by  ordinary  particle  scattering. 

Although  the  use  of  fluorescence  rather  than  scattering  has  many 
obvious  advantages,  particularly  for  backscatter  systems,  no  reports 
of  the  use  of  this  technique  have  appeared  in  the  literature  since 
the  original  publication  of  Stevenson  et  al.  in  1975.  Private  com¬ 
munication  with  Dr.  T.  Powell  of  NASA-Lewis  and  Dr.  Carl  Williams  of 
United  Technology  Research  Laboratories  have  confirmed  that  the  fluores¬ 
cence  method  exhibits  the  advantages  cited  previously.  However,  they  have 
encountered  one  practical  difficulty  with  the  method  described  by  Stevenson 
et  al.,  a  difficulty  that  our  system  will  avoid.- 

Stevenson  et  al.  (1975)  produced  fluorescent  particles  by  atomizing 
liquids  containing  appropriate  fluorescent  dyes.  These  liquid  droplets 
are  suitable  for  air  flows  but  are  unsuitable  for  measurements  in  liquids. 
Even  in  gas  streams,  the  droplets  tend  to  impinge  on  surfaces,  leaving  a 
fluorescent  '-esidue  on  the  surface  that  produces  background  signals  similar 
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to  those  caused  by  wall  imperfections  in  conventional  scattering  systems. 
According  to  Williams  (private  communication) ,  this  prohibits  the  use  of 
the  fluorescence  method  within  3-4  probe  volumes  of  the  surface. 

Our  research  will  employ  solid  particles  containing  fluorescent  dye, 
thereby  eliminating  entirely  the  difficulties  encountered  with  the  liquid 
droplet  system  used  by  Stevenson  and  others. 

In  summary,  a  backscatter  LDA  system  using  fringe-mode  single¬ 
particle  detection  is  found  to  be  superior  to  other  optical  arrangements 
for  velocity  measurements  near  solid  surfaces.  The  use  of  solid  fluores¬ 
cent  particles  will  eliminate  the  large  background  signals  that  have 
characterized  previous  experiments  using  scattered  laser  light  and  liquid 
fluorescent  particles. 

B .  System  Design 

The  operation  of  the  laser  Doppler  velocimeter  in  the  "fringe  mode" 
is  based  on  fairly  simple  concepts.  If  two  coherent,  linearly  polarized 
laser  beams  are  brought  to  a  focus  at  a  common  point,  a  set  of  intensity 
fringes  is  produced  at  the  beam  intersection.  These  fringes  are  equally 
distant  from  each  other;  their  spacing  can  be  calculated  from  the  geometry 
of  the  optical  train.  A  particle  traversing  this  field  of  spatially 
varying  light  intensity  scatters  light  with  an  intensity  proportional  to 
that  of  the  incident  laser  intensity.  The  scattered  light  is  collected 
by  a  photodetector  and  the  signal  is  transmitted  to  a  counter.  The  counter 
clocks  the  time  taken  by  the  particle  to  cross  the  fringes  and  simultaneously 
determines  the  frequency  of  crossing.  The  particle  velocity  perpendicular 
to  the  fringe  planes  is  calculated  by  taking  the  product  of  the  frequency 
and  the  fringe  spacing. 

In  order  to  complete  a  quantitative  analysis  of  the  laser  scattering 
process,  it  is  necessary  to  describe  at  some  length  the  fringe  field  pro¬ 
duced  by  two  focused,  intersecting,  coherent  laser  beams.  The  following 
sections  explain  the  interference  phenomenon,  the  nature  of  Gaussian  beams, 
and  the  details  of  the  region  produced  by  the  intersection  of  the  laser 
beams.  Tracer  particles  are  then  selected  for  their  abilities  to  follow 
the  flow  and  produce  a  satisfactory  signal.  The  concluding  section  of 


tnis  chapter  describes  the  shot  noise  associated  with  background  scatter. 
Interference  and  Fringe  Spacing 

When  two  light  beams  originating  from  the  same  source  are  super¬ 
imposed,  interference  results.  This  phenomenon  of  interference  is 
characterized  by  intensity  variations  where  the  maxima  may  exceed  the 
linear  sum  of  the  intensities  of  the  two  beams  and  where  the  minima  may 
be  zero.  In  the  case  of  plane,  monochromatic,  linearly-polarized  light 
waves,  this  interference  is  seen  as  a  series  of  bright  and  dark  fringes 
which  represent  regions  of  maximum  and  minimum  intensity,  respecti vely. 

The  output  of  a  laser  can  be  ideally  described  as  a  plane,  mono¬ 
chromatic,  linearly-polarized  light  beam.  In  LDA,  the  output  beam  is 
split  by  means  of  an  optical  beamsplitter  into  two  beams,  each  with  a 
lower  intensity  than  the  original  beam.  The  two  resulting  beams  are  sub¬ 
sequently  intersected  with  one  another  to  form  a  region  of  interference 
that  is  characterized  by  a  series  of  bright  and  dark  intensity  fringes. 

Figure  9  depicts  two  such  beams,  represented  by  electric  field 
vectors  e-j  and  e 2,  intersecting.  ^  an<J  n2  repreSent  the  directions  of 
propagation,  and  the  lines  orthogonal  to  these  directions  represent  the 
locations  of  a  crest  or  a  peak  in  the  wave,  and  can  be  expressed 
as  complex,  scalar  quantities;  the  scalar  representation  arises  because 
both  beams  are  not  only  plane  and  monochromati c,  but  linearly  polarized. 

el  =  (Eo}l  exp  -  kh 2  +  <t>])3  (43a) 


e2  =  (Eo^2  exp  ^'(27rvt  "  +  *2>1  t43b* 

In  the  above  equations,  (E0).j  represents  the  amplitude  of  the  elec¬ 
tric  field  vector  for  the  i—  beam,  <  represents  the  wave  number  and 
equals  2tt/x,  and  represents  the  phase  of  the  i—  beam.  The  total 
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electric  field  in  the  intersection  region  can  be  expressed  as: 

c  -  +  c2  (44) 

The  intensity  variation  across  the  region  is  given  by  the  time 
average  of  tne  square  of  the  electric  field,  where  the  time  average 
is  taken  over  a  large  interval  of  time  compared  with  the  period,  2ir/v. 
The  following  expression  is  obtained  for  intensity  (Durst,  Melling, 
and  Whitelaw,  1976): 

1  =  +  <Eo}2  +  C0S  "  r'2)  +  '  *2)]  <45) 

Assuming  that  the  phases  of  the  beams  are  equal  to  each  other  at  t  =  0, 
the  maximum  intensity  occurs  when  (c2  -  n2)  *  Nx,  or  at  an  integral 
number  of  wavelengths.  The  minimum  intensity  occurs  when  (?2  -  n2)  * 

(N  +  1)X. 

Figure  10  depicts  an  enlarged  view  of  the  interference  region 
where  the  dark  horizontal  lines  represent  the  intersections  of  the 
crests  in  the  beams.  The  fringe  spacing  can  be  calculated  with  the 
aid  of  this  diagram  and  a  coordinate  transformation.  From  the  geometry 
in  Figure  10. 

n 2  =  x2  cos  3  ~  xi  S1n  9 

s2  =  x2  cos  e  +  x.j  sin  e 
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where  9  is  tr.e  half-angle  of  intersection.  Subtracting  n2  ^roir'  c2 


gives: 


(c2  -  n2)  =  2x1  sin  9 


the  fringe  is  given  by: 

U2  -  n2)  =  NX  =  2(x-j )N  sin  e 


and  the  (N  +  1 )—  fringe  by: 

(?2  -  n2)  «  (N  +  1)X  =  2(x1)n+1  sin  9 


Hence, 


/y  \  _  wa  .  /x  \  = 

1 Vn  2  sin  e  ’  v  rN+1  2  sin  0 

and  the  distance  between  the  two  fringes,  d^,  is: 

df  =  ^-Pn+I  "  ^X1^N  ”  2  sin  6 

Thus,  the  fringe  separation  is  dependent  only  on  the  wavelength  of  the 
intersecting  beams  and  the  angle  at  which  they  intersect. 

Rudd's  Fringe  Model 

The  principle  of  the  interference  fringe  model  proposed  by  Rudd  in 
1961  [see  Durst  et^  al_.  (1976)]  provides  a  clear  physical  insight  into  the 
physics  of  LDA  and  a  quantitative  basis  for  analyzing  laser  Doppler  signals. 
A  particle  passing  through  a  set  of  interference  fringes,  as  seen  earlier 
in  Figure  10,  scatters  light  with  an  intensity  proportional  to  that  of  the 
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light  incident  upon  it.  The  scattered  light  is  detected  by  a  photo¬ 
multiplier  which  converts  the  incoming  optical  signal  into  an  electrical 
signal.  The  electrical  signal  is  passed  on  to  an  electrical  counter  which 
calculates  the  frequency  of  crossing,  This  signal  frequency  is  re¬ 
lated  to  the  velocity  of  the  particle  in  the  direction  perpendicular  to 
the  fringe  planes,  Vx,  by  the  following  equality: 


Vx  2VX  sin  e 
VD  "  37  - X - - 


(48) 


where  was  calculated  in  Equation  (47). 

Vp  is  commonly  referred  to  as  the  Doppler  frequency,  although  when 
speaking  in  terms  of  Rudd's  fringe  model  the  term  Doppler  frequency  is  not 
quite  accurate.  The  above  expression  for  however,  is  the  same  as  that 
obtained  for  the  signal  frequency  based  on  Doppler  considerations  (Durst, 
Melling,  and  Whitelaw,  1976).  The  semantics  problem  of  whether  or  not  the 
fringe  system  just  described  can  justifiably  be  called  a  Doppler  system  is 
not  of  practical  interest.  Further  reference  to  the  Doppler  anemometer  and 
the  Doppler  frequency,  vQ,  will  be  understood  to  comprise  the  above  mentioned 
fringe  system  and  signal  frequency,  respectively. 

Gaussian  Beams 

In  LDA  systems  the  laser  is  almost  always  operated  in  the  TEMQ0 
(fundamental  transverse  electromagnetic)  mode  resulting  in  a  beam  with 
a  Gaussian  intensity  profile.  Figure  11  depicts  such  a  Gaussian  beam 
with  a  beam  waist  of  2wQ,  a  beam  radius  of  w(z)  at  some  arbitrary  z-location, 
and  plane  wavefronts  with  radius  of  curvature,  R(z).  The  beam  waist,  which 

is  the  narrowest  portion  of  the  beam,  and  the  beam  radius,  w(z),  are  defined 

2 

at  the  1/e  intensity  points.  w(z)  is  given  by: 


w2U)  *  w02  [1  +  (A*-)2] 


(49) 
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Tne  radius  of  curvature,  R(z),  of  the  wavefront  can  be  computed: 


R(z)  =  z  [1  +  { 


ttw  2  2 

— )  3 

xz  ‘  J 


(50) 


From  this  equation  it  is  seen  that  at  z  *  0,  R(z)  goes  to  infinity 
and  the  wavefronts  are  essentially  planar. 

The  original  size  and  location  of  the  beam  waist  is  determined 
by  the  configuration  of  the  laser  cavity.  The  size  and  location  of 
the  beam  waist  after  being  focused  by  a  lens  (see  Figure  12)  are  given 
by: 


Z1  =  f  + 


(zo  ~  f)2  *  (1TW°2)2 


(51) 


_ wof 

I  nw  2  2  1/2 

[<*o  -  f>  +  (-T-J  3 


(52) 


where  f  is  the  focal  length  of  the  lens.  Hence  the  focused  beam  waist  is 
generally  not  located  in  the  front  focal  plane  of  the  lens.  The  exceptional 
case  occurs  when  zQ  is  in  the  back  focal  plane. 

When  two  Gaussian  beams  are  focused  by  a  lens  of  focal  length  f,  the 
volume  circumscribed  by  the  1/e2  intensity  points  of  the  fringes  is  called 
the  probe  volume.  Since  the  beam  waists  are  not  found  at  the  front  focal 
plane  of  the  lens  (except  when  zQ  =  f)  they  do  not  intersect  at  the  center 
of  the  probe  volume  and  may  even  be  outside  the  probe  volume.  In  either  case 
some,  if  not  all,  of  the  wavefronts  within  the  probe  volume  will  have  a 
finite  radius  of  curvature  and  will  exhibit  a  varying  fringe  spacing  as 
shown  in  Figure  13.  This  results  in  signals  of  varying  frequencies  depending 
on  the  path  traversed  by  the  particle  through  the  probe  volume.  Durst  and 


-55- 


Stevenson  (1979)  found  the  gradient  of  the  Doppler  frequency  along  the 
probe  volume  to  be 


(53) 


To  assure  that  the  waists  of  both  beams  intersect,  an  additional 
optical  component  consisting  of  both  a  positive  and  a  negative  lens  is 
positioned  between  the  laser  and  the  transmitting  lens.  The  position 
of  the  beam  waists  relative  to  the  lens  can  then  be  varied  and  made  to 
occur  at  the  front  focal  plane  of  the  lens. 

Calculation  of  Probe  Volume  Dimensions 

If  the  beam  waists  do  intersect  after  being  focused  by  a  lens,  an 
elliptical  probe  volume  is  formed.  The  dimensions  of  this  probe  volume 
are  determined  by  the  geometry  of  the  optical  system. 

Figures  14  and  15  depict  a  fairly  complicated  set-up  consisting  of 
two  beams  expanded  in  air  and  subsequently  focused  by  two  lenses  into  a 
medium  with  an  index  of  refraction  different  from  that  of  air.  Figure  14 
is  a  side  view  of  the  set-up;  consequently  only  one  of  the  beams  is  seen. 
Figure  15  illustrates  the  position  of  the  other  beam  relative  to  the  first. 
The  result  is  an  ellipsoid  consisting  of  a  set  of  parallel  fringe  planes 
equally  distant  from  one  another.  The  diameter,  dm,  length,  l  >  and  width, 
wm,  can  be  calculated  given  i,  r  ,  dbQ,  ex,  f j ,  fj,  l,  and  z,  where: 


i  =  index  of  refraction  of  medium 

rQ  =  distance  between  beams  and  optical  axis 

dbo  =  diameter  of  beams  entering  beam  expander 
(assumed  to  be  diameter  leaving  laser) 

ex  *  beam  expansion  ratio 


Ai 
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fj,  f2  3  focal  lengths  of  lenses  1  and  2,  respectively 

i,  z,  =  variable  distances  representing  locations  of 
1  lens  2  and  medium's  boundary  with  respect  to 
lens  1 

In  Figure  14,  the  beams  enter  the  beam  expander  with  an  expansion 
ratio  of  ex,  at  a  distance  of  r  ,  and  a  diameter  of  <V  They  leave  at 
a  distance  of  r^ ,  and  a  diameter  of  dbl  where: 

rl  ■  ex  ro 
dM  *  ex  "to 

If  several  beam  expanders  are  used  in  series,  the  total  beam  expansion 
ratio  is  the  product  of  the  individual  expansion  ratios. 

6,,  the  angle  which  would  exist  between  the  beams  and  the  optical 
axis  if  lens  2  and  the  medium  were  absent  is  given  by: 

-1  rl 

*  tan 

The  beams  enter  lens  2  at  a  distance  of  r 2  from  the  optical 
axis,  defined  by: 

r2  *  £(tan  e^)  . 

feff’  the  effective  focal  length  with  respect  to  lenses  1  and  2,  and  e2, 
the  angle  between  the  beams  and  the  optical  axis,  are  found  by  using  the 
lens  equation  and  noting  the  geometry: 


Figure  15.  Isometric  View  of  Beam  Paths  for  Measurement  of  One  Velocity 
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f  ,f2  ' 

feff  (fg  +  z) 

-1  r? 

92  *  tan  {T—) 
eff 

Using  Snell's  law  of  refraction  and  assuming  the  index  of  refraction 

l 

of  air  to  equal  one,  the  angle,  82,  at  which  the  beams  enter  the  medium, 
and  X‘ ,  the  wavelength  at  which  they  propagate  through  the  medium  are  cal¬ 
culated  to  be: 


I 


,  sin  6, 
sin'1  (-— l 


V 


\ 

T 

1 


i.',  the  location  of  the  probe  volume  before  refraction  by  the 
medium  is  given  by: 


*'  a  Z1  '  "  feff) 

This  gives  rise  to  the  determination  of  f  ,  the  effective  focal  length 
with  respect  to  lens  1,  lens  2,  and  the  medium,  and  db2,  the  diameter  of 
the  beam  entering  the  medium: 


tan  e2‘ 


1 1  tan  e 
tan  e2‘ 


j  1  £'  ,  £' 

“  ■  feff  db2  * 


From  Figure  16,  A',  the  angle  of  the  "cone"  formed  by  the  focused 
beam,  is  found  (assuming  small  A'): 


* 
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A* 


COS  9  2 

- r— 

eff 


I 


The  expression  for  d-j  is  given  by: 


d 


1 


=  4  XV 
IT  A' 


<2  ,  the  half  angle  between  beams  1  and  2,  and  beams  1  and  3,  is  found 
through  geometry: 


K 


2 


sin 


-1 


sin  e  ' 

( - -J-) 


and  the  probe  volume  dimensions  as  seen  in  Figure  16  are: 


d 


m 


COS  Kg' 


m 


sin  Kg' 


The  fringes  within  the  probe  volume  are  plane,  parallel,  and  equi 
distant.  The  fringe  spacing  was  calculated  earlier  in  Equation  (47)  as 


d 


f 


x 

2  sin  Kg1 


Hence,  the  number  of  fringes,  N^,  in  the  probe  volume  can  easily  be 
calculated: 


Figure  17  depicts  the  LOA  optical  train  to  be  used  to  measure  one 
velocity  component.  The  beam  from  the  laser  is  split  in  two  and  subse¬ 
quently  focused  by  a  converging  lens  into  an  ellipsoidal  probe  volume 
with  parallel  fringe  planes.  The  Doppler  signals  from  particles  tra¬ 
versing  the  probe  volume  are  collected  by  the  photomultiplier  and  passed 
to  the  electronic  counter.  The  analog  output  from  the  photomultiplier 
is  also  viewed  using  an  oscilloscope.  The  digital  output  is  sent  to  the 
computer,  where  it  is  stored,  via  an  interface  designed  and  built  in  this 
laboratory  by  G.  Muntz,  D.  Dewey,  and  C.F.  Dewey,  Jr. 

All  of  the  optical  components,  except  for  that  part  of  the  retarder 
attached  to  the  laser,  are  physically  coupled  to  one  another  creating  an 
optical  train. 

The  various  components  of  the  transmitter  portion  of  the  optical 
train  shown  in  Figure  17  were  chosen  to  produce  a  variety  of  probe  volume 
dimensions  and  fringe  spacings  within  the  probe  volumes.  With  two  differ¬ 
ent  front  lenses  with  focal  lengths,  f,  of  160  mm  and  80  irm,  and  using 
prism  spacings,  VQ,  (the  inner  prisms  of  the  beam  translator)  of  6.5  mm 
to  19.5  mm  from  the  optical  axis,  the  operating  conditions  given  in  Table  2 
can  be  achieved. 


LDA  Set-up  for  Backscatter  Mode,  One  Velocity  Component. 
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TABLE  2.  PROBE  VOLUME  DIMENSIONS1*3  (SIMP1  E  SET-UP2! 

See  Fig.  16  for  symbol  definitions. _ 

~Z  FT  ~  TZ  31  vT  dL  n7 
t  o  ex  m  m  m  ff 

(mm)  (mm)  {ym)  (ym)  (ym)  (ym) 


160 

6.5 

1 .0 

2555.74 

73.39 

73.36 

8.96 

8 

150 

6.5 

1.9 

711.01 

38.75 

38.69 

4.73 

8 

160 

13.0 

1.0 

1284.19 

73.66 

73.54 

4.49 

16 

160 

13.0 

1.9 

361.82 

39.26 

39.03 

2.38 

16 

150 

19.5 

1.0 

863.14 

74.11 

73.84 

3.01 

24 

30 

6.5 

1.0 

642.09 

36.83 

36.77 

4.49 

8 

80 

6.5 

1.9 

180.91 

19.63 

19.52 

2.38 

8 

80 

13.0 

1.0 

327.36 

37.37 

37.13 

2.27 

16 

80 

13.0 

1.9 

96.77 

20.64 

20.19 

1.23 

16 

80 

19.5 

1.0 

225.26 

38.26 

37.72 

1.54 

24 

^Calculated  for  a  beam  diameter  of  1.43  rm  leaving  the  beam 
waist  adjuster  and  a  medium  with  index  of  refraction  =1. 

p 

A  more  complicated  set-up  will  be  examined  in  later  calculations. 

3 

ex  is  the  expansion  ratio  of  the  beam  diameter 

In  Figure  17  the  transmitted  beams  are  indicated  by  the  dark  solid 
lines.  The  paths  go  through  a  retarder,  a  beam  waist  adjuster,  a  neutral 
beamsplitter,  a  beam  displacer,  a  beam  translator,  a  beam  expander,  and  a 
front  lens.  The  retarder  consists  of  two  auarter-wave  plates;  the  first 
plate  circularly  polarizes  the  incoming  laser  beam  and  the  second  linearly 
polarizes  it  in  the  correct  polarization  direction  relative  to  the  optical 
axis.  The  polarization  rotates  with  the  optical  train  due  to  the  physical 
attachment  between  the  second  quarter-wave  plate  and  the  afore  mentioned 
train. 

The  beam  waist  adjuster  causes  the  beam  waist  of  the  two  focused  beams 
to  be  located  at  the  focal  point  of  the  front  lens.  The  neutral  beamsplitter 
splits  the  incoming  beam  into  a  50:50  intensity  ratio.  The  beam  displacer 
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shifts  the  beam  lying  on  the  optical  axis  a  fixed  radial  distance  away. 

When  this  is  done,  the  radii  from  the  optical  axis  to  the  new  beam  loca¬ 
tions  form  right  angles  with  each  other.  The  beam  translator  allows  the 
length  of  the  radii  to  the  beams  to  be  varied  and  consequently  affects 
both  the  diemsnions  of  the  probe  volume  and  the  number  and  spacing  of 
the  fringes.  The  beam  expander  is  used  in  cases  where  it  is  desirable 
to  decrease  tne  dimensions  of  the  probe  volume  by  increasing  the  inter¬ 
section  angle  of  the  beams.  The  final  component  in  the  transmitter  path, 
the  front  lens,  is  of  prime  importance  since  it  not  only  focuses  the  beams 

to  form  a  probe  volume  but  acts  as  a  collection  agent  in  the  receiving  path. 

The  lenses  must  be  as  free  from  aberrations  as  possible  to  ensure  high 
quality  signals.  In  addition,  it  is  advantageous  to  employ  achromatic 
lenses  to  correct  for  aberrations  arising  from  the  variation  of  the  refrac¬ 
tive  index  of  glass  with  wavelength.  This  is  penally  important  for  work 
with  fluorescent  signals  where  the  excitation  wavelength  differs  from  the 
signal  wavelength. 

The  Optical  Receiver 

The  receiver  portion  of  the  optical  train  collects  the  backscattered 

(or  emitted)  light  from  the  particles  and  hence  contains  the  sought  after 

velocity  information.  It  consists  of  the  front  lens,  the  beam  expander, 
the  pinhole  section,  the  backscatter  section,  the  photomultiplier  optics, 
and  the  photomultiplier  tube.  The  receiver  path  is  denoted  in  Figure  17 
by  the  shaded  area. 

The  front  lens  collects  the  backscattered  or  emitted  light  directly, 
making  it  necessary  to  have  it  coated  with  an  anti -reflection  film.  After 
collection  by  the  front  lens,  the  light  is  passed  through  the  beam  expander. 
The  beam  expander  section  is  used  in  conjunction  with  the  beam  translator  to 
increase  the  convergence  angle  of  the  transmitted  beams.  It  does  not  alter 
the  collection  efficiency  or  effective  collection  angle  of  the  optical  re¬ 
ceiver. 

The  front  lens  of  the  pinhole  section  focuses  the  collection  rays 
through  a  tiny  pinhole  removing  a  large  fraction  of  the  unwanted  light 
generated  through  spurious  reflections.  The  rays  are  sent  to  the  backscatter 
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section  consisting  of  a  plane  mirror  angled  at  45°  to  the  optical  axis. 

The  mirror  directs  the  light  through  the  photomultiplier  optics  where 
they  are  focused  through  another  pinhole.  This  pinhole  acts  as  a  spatial 
filter  further  eliminating  light  scattered  or  reflected  from  the  various 
optical  surfaces. 

The  light  is  finally  conveyed  to  the  photomultiplier  where  the 
incoming  optical  signal  is  converted  to  an  electrical  signal. 

Versatility  of  Optical  Components 

After  analyzing  the  performance  and  prices  of  the  two  major  com¬ 
panies  specializing  in  IDA  apparatus,  the  optics  and  electronics  were 
ordered  from  DISA  Electronics.  Careful  consideration  was  given  to  both 
the  quality  of  the  components  and  the  flexibility  of  the  entire  system. 

Since  the  immediate  objective  of  the  experimental  work  in  this  project 
was  to  determine  whether  or  not  the  use  of  fluorescent  particles  as 
tracers  worked  in  an  LDA  set-up,  the  initial  system  had  to  be  kept  simple. 

The  ultimate  objective,  however,  is  to  be  able  to  measure  shear  stresses  in 
two-dimensional  flows,  perhaps  with  turbulent  characteristics. 

In  such  a  case  it  is  necessary  to  create  an  orthogonal  fringe  system 
in  as  small  a  probe  volume  as  possible,  with  frequency  shifting  of  one  or 
more  beams  to  distinguish  a  positive  velocity  from  a  negative.  With  this 
in  mind,  the  set-up  depicted  in  Figure  18  was  envisioned.  In  this  set-up, 
a  neutral  beamsplitter  is  first  employed  to  split  the  laser  beam  (cyan- 
colored  if  a  two-color  set-up  is  desired)  into  equal  intensity  ratios. 

In  the  mono-color  set-up,  another  neutral  beamsplitter  fol 'ows  the  first; 
n  the  two-color  set-up,  the  neutral  beamsplitter  is  followed  by  a  dichroic 
beamsplitter.  The  second  beamsplitter  divides  only  one  of  the  beams,  either 
into  equal  intensity  ratios  in  the  mono-color  case  or  into  a  green  beam  and 
a  blue  beam  in  the  two-color  case. 

With  the  two-color  set-up,  the  cyan  and  the  blue  beams  intersect  to 
form  a  series  of  parallel  fringe  planes  orthogonal  to  the  fringe  planes 
formed  by  the  cyan  and  the  green  beams.  Similarly,  in  the  mono-color  set-up, 
the  most  intense  beam  combines  with  each  of  the  two  beams  to  form  two  ortho¬ 
gonal  sets  of  fringe  planes.  A  two-color  backscatter  system  can  distinguish 
between  scattering  from  the  two  orthogonal  fringe  sets  on  the  basis  of  wave 
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length.  With  the  mono-color  system  (i.e.,  a  single  wavelength),  Bragg 
cells  must  be  used  to  shift  the  apparent  laser  frequency  to  obtain  dis¬ 
crimination  between  the  two  velocity  dimensions. 

Bragg  cells  are  usually  employed  to  optically  shift  one  beam  at 
a  time  by  40  MHz.  This  produces  a  set  or  sets  of  moving  fringe  planes. 

A  stationary  particle  situated  in  the  probe  volume  produces  a  frequency 
equal  to  the  shift  frequency.  Movement  in  the  same  direction  as  the 
fringes  lowers  the  signal  frequency  and  in  the  opposite  direction  raises 
it.  Both  particle  speed  and  direction  can  subsequently  be  deduced. 

Two  beam  expanders  in  combination  with  lenses  of  high  f-number 
(focal  length/diameter)  significantly  reduce  the  dimensions  of  the  probe 
volume.  Table  3  displays  the  dimensions  of  the  probe  volume  as  well  as 
the  number  and  spacing  of  the  fringes  following  the  calculations  described 
previously. 


TABLE  3.  PROBE  VOLUME  DIMENSIONS1 

(COMPLICATED  SET-UP2! 

See  Fig  16. 

r 

o  ex 

(mm) 

- 1 - 

m 

(urn) 

dm 

(urn) 

w  d,  N£ 

m  T 

(um)  (urn) 

6.5 

1.0 

548.47 

29.51 

29.46 

3.60 

8 

6.5 

1.9 

156.07 

15.80 

15.72 

1.92 

8 

6.5 

1.92 

47.27 

15.99 

8.65 

1 .06 

15 

13.0 

1.0 

282.63 

30.08 

29.91 

1 .83 

16 

13.0 

1.9 

86.44 

16.86 

16.55 

1 .01 

16 

19.5 

1.0 

197.76 

31.02 

30.64 

1 .25 

24 

Calculated  for  a  beam  diameter  of  1.43  mm  leaving  the 
beam  waist  adjuster  and  a  medium  with  index  of  refrac 
tion  =  1.33. 

2 

Front  lenses  with  f,  =  160  mm,  f?  *  80  mm;  Distance, 
l  *  120  mm;  Distance,  Zp  *  95  mnr(see  Figure  14). 


Geometrical  Configuration 

Figure  19  is  a  schematic  of  the  geometrical  configuration  of  the 
laser  and  optics.  Both  laser  and  optics  had  to  be  securely  fastened  to 
a  vi brat ion -free  base  since  the  smallest  movement  resulted  in  misalign¬ 
ment  and  a  distortion  in  or  even  loss  of  the  fringe  pattern  in  the  probe 
volume.  They  were  bolted  to  a  2-1/2'  x  4*  granitetable  sitting  on  shock- 
absorbent  pads.  The  pads,  purchased  from  Air-Loc^,  were  designed  to 
withstand  loads  up  to  1500  lbs  and  absorb  85%  of  the  vibration. 

The  mirrors  were  purchased  from  Newport  Research  Corporation  and 
were  rated  99‘i  reflective  at  visible  wavelengths.  The  optical  train 
rested  on  an  optical  rail  to  permit  ease  in  translation  of  the  components 

C.  Continuing  Research 


Fabrication  of  the  optical  stands  and  integration  of  the  optical 
system  are  the  next  objectives.  Following  this,  measurements  will  be 
made  with  targets  of  known  velocity  to  calibrate  the  LDA  system.  Then 
measurements  wi’l  be  made  using  fluorescent  particles  and  an  estimate  of 
signal-to-noise  ratio  will  be  made.  Finally,  the  LDA  will  be  applied  to 
practical  shear  flows  near  solid  boundaris  where,  heretofore,  LDA  has  not 
been  a  practical  metnod  of  measurement. 
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IV.  SUMMARY 


Fluid  shear  stress  is  one  of  the  most  important  quantities 
characterizing  the  interaction  between  a  flowing  fluid  and  a  solid  body. 

A  knowledge  of  surface  shear  stress  is  essential  in  determining  the  drag 
on  submerged  bodies,  the  energy  losses  in  internal  flows,  local  surface 
heat  transfer  rates  in  both  internal  and  external  flows,  flow  separation, 
cavitation,  and  unsteady  flow  behavior  in  rotating  machinery. 

In  spite  of  many  previous  investigations  of  both  laser  Doppler 
anemometry  (LDA)  and  surface  hot-film  gauges,  the  use  of  these  methods 
for  determining  wall  shear  stress  is  poorly  developed.  Conventional 
LDA  measurements  near  surfaces  are  plagued  by  spurious  background  scattering 
and  marginal  resolution.  Use  of  surface  hot-film  gauges  is  severely  in¬ 
hibited  by  the  necessity  of  empirical  calibration  and  an  incomplete  under¬ 
standing  of  their  response  in  turbulent  and  unsteady  flows. 

This  research  program  will  develop  new  methodology  for  the  use  of 
these  two  complementary  techniques  as  quantitative  measurement  tools  for 
determining  local  surface  shear  stress.  The  LDA  system  will  be  based  on 
detecting  laser-induced  fluorescence  from  solid  particles  in  a  so-called 
"backscatter"  mode.  Two  individual  velocity  components  will  be  measured 
by  using  multiplexing  in  the  frequency  domain.  The  hot-film  probes  will 
employ  multiple  elements  such  that  conventional  calibration  procedures 
can  be  eliminated.  Both  theory  and  experiment  are  being  used  to  determine 
the  utility  of  LDA  and  hot-film  instrumentation  in  unsteady  flows. 

Our  analysis  of  hot-film  probe  performance  considers  an  arbitrarily 
shaped  probe  (or  set  of  probes)  flush-mounted  on  a  plane  thermally  conduc¬ 
ting  substrate,  possessing  an  arbitrary  heat  flux  distribution,  and  exposed 
to  a  steady  uniform  fluid  shear.  We  have  formulated  the  complete,  three- 
dimensional  fluid  and  substrate  energy  equations  and  have  coupled  them 
through  their  associated  boundary  conditions.  A  novel  solution  procedure 
that  combines  Fourier  transformation  of  the  governing  equations,  their 
analytical  solution  in  transform  space,  and,  finally,  numerical  inversion 
to  predict  physical  space  temperature  distributions  and  heat  fluxes  is  being 


developed.  This  approach  has  some  clear  potential  advantaaes  --  as  well 
as  some  possbile  weaknesses  --  that  are  being  systematically  explored. 

One  important  advantage  is  that  the  solutions  obtained  can  be  limited  to 
any  plane  of  interest  (particularly  the  fluid-solid  interface);  solutions 
for  temperature  and  heat  flux  distributions  elsewhere  in  the  fluid  or  sub¬ 
strate  need  not  be  generated  unless  desired.  This  method  should  thus  prove 
to  be  computationally  efficient  for  determining  what  is  of  greatest  interest 
in  probe  analysis  -  the  probe  temperature  and  heat  flux  distributions.  Since 
the  full  three-dimensional  fluid  and  solid  energy  equations  are  considered, 
the  technique  should  be  able  to  handle  a  wide  range  of  fluids  and  flow 
speeds,  probe  geometries  and  substrate  properties.  The  use  of  this  tech¬ 
nique  will  enable  us  to  obtain  detailed  descriptions  of  the  temperature  fields 
and  heat  flows  associated  with  a  given  probe  design.  It  will  also  enable 
us  to  investigate  the  sensitivity  of  probe  performance  to  various  flow  and 
substrate  parameters  and  to  thus  determine  a  set  of  "optimum"  designs  which 
could  then  be  "numerically  calibrated".  Potential  weaknesses  of  this  appraoch 
approach  may  be  associated  with  numerical  difficulties  in  obtaining  highly 
accurate,  yet  efficient,  solutions  to  the  resulting  two-dimensional  Fourier 
inversion  integrals  for  desired  probe  configurations. 

Our  laser  Doppler  anemometry  (LDA)  efforts  are  keyed  to  the  use  of 
fluorescent  tracer  particles  in  combination  with  optical  filters  to  overcome 
the  severe  limitations  of  conventional  LDA  techniques  near  solid  surfaces. 
Here,  background  scattering  is  large  and  LDA  becomes  very  difficult  because 
the  "shot  noise"  (statistical  uncertainty)  from  the  background  swamps  the 
signal  arising  from  tracer  particles  moving  through  the  measurement  volume. 

The  fluorescent  emission  is  at  longer  wavelengths  than  the  incident 
absorbed  wavelength;  consequently,  optical  filters  can  be  employed  to  filter 
out  the  background  scattering  while  passing  the  fluorescent  signals.  The 
theoretical  estimates  we  have  made  predict  large  (factor  of  1000)  increases 
in  the  signal-to-noise  ratio  using  the  fluorescent  tracer  particles  over  that 
using  conventional  scattering  particles. 
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These  methods  may  be  used  with  conventional  backscatter  optics 
and  either  one  or  two  component  systems;  however,  multiplexing  in  the 
Bragg  shift  frequency  domain  rather  than  in  the  optical  (wavelength) 
domain  is  required  because  the  fluorescence  is  at  a  wavelength  longer 
than  the  excitation  wavelength. 
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